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Preface 


The primary objective of this Interpretation Techniques Development 
Study was to quantitatively relate the photoqraphfc characteristics of the 
ERTS-I images to measurements of reflectance as well as geologic terrain 
classifications, as exist on geologic maps, in south eastern Arizona, 

Ground spectral reflectance measurc'ment;;< of the V«llScox Playa were 
made on two different dates* These data collected at twelve stations within 
the Playa were made on the same day as ERTS“1 coverage. Simultaneous soil 
samples analyzed for moisture content showed that change in spectral ref- 
lectance was not correlated with soil moisture but rather represents standing 
surface water. 

Reflectance measurements of three tree species in the Prescott National 
Forestwere collected prior to the launch of ERTS-I, These reflectance 
spectra could not be correlated with ERTS-1 images at locations where known 
plantations of These trees existed. 

A computerized "System 800" was used to quantitatively examine ERTS-I 
image density characteristics. This system used TV digitization of four 
ERTS-I photographic images in registration into 294,912 pixels for each photo, 
each pixel having 8 bits (256 levels). The data were stored in a computer 
disc, and processed by PDP-11 computer using a variety of statistical and 
display programs. Output was on a 8 color TV display. This system was used 
to make both spectra I and temporal analysis of the Willcox Playa as well as 
the PrescotT National Forest and Phoenix, Arizona, Relative atmospheric 
effects in each of the four bands weredetermi ned along with the statistical 
characteristics of the brightness In each ERTS-1 band. These data were used 
to develop procedures to correct for atmospheric effects in each of the 
four bands. 

Photo-optical additive color techniques v<ere used to perform geologic 
interpretations using ERTS-1 black-and-white positives. Processing pro- 
cedures were developed which when used to make positives from standard 
release negatives, when these positives were used In an additive color 
viewer, produced optimum color Image characteristics of soils and rocks. 

Many features appearing in the imagery were found not to appear on the 
most recent geologic maps. 

Additive color viewer procedures were developed which produced tem- 
poral changes of standing surface water between different ERTS-1 passes. 

These procedures were adopted by the U. S, Geological Survey, Water Resources 
Division as a standard procedure to be used in mapping floods using ERTS, 

The 1973 Mississippi Flood was mapped using the temporal additive color 
techniques developed in this study. 
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Section i 


Introduction 

This report presents the results of nn Interpretation Techniques 
Development study, the initial objective of which was to quantitatively 
relate photographic characteristics of the ERTS-I images to measure- 
ments of the reflectance, and of the characteristics of ground objects. 
Positive and negative Multi spectral Scanner (MSS) images taken prin- 
cipally over the Willcox Playa, Arizona as received from Goddard Space 
Flight Center were the basic input to this project. 

Initial Data Analysis 

An analysis of the ERTS-I photographic characteristics was per- 
formed initially, the purpose being to establish optimum procedures 
for printing positives of semi-arid desert-like terrain. A procedure 
was developed which when used to print positive transparencies from 
NASA ERTS-I release negatives produced excellent composite color images 
when projected on an Additive Color Viewer screen. 

An extensive qualitative interpretative analysis was performed 
using an Additive Color Viewer and the reprocessed positive images* 

The purpose here was to examine spectral reflectance image character- 
istics of a large and relatively homogenous object such as the 
Wi 1 1 cox Pl aya . 

Ground Truth 

The principal area of interest In this study, the Willcox Playa, 

Is a large dry lake located in southeastern Arizona* The central part 
of the Playa has a hard and compact crust* The outer areas of the 

r 


PAG® ® 




Playa contain soft, dry and porous surfaces v/hich become flooded durinq 
the infrequent rains. 

Reflectance spectra were taken in situ on the Playa before the 
launch of ERTS-I and during 1973, These latter measurements were made 
at twelve sample stations in the Playa and correlated with soil moisture 
measurements simultaneously taken. 

Reflectance measurements were also taken at the Prescott National 
Forest In north central Arizona. Three plant species were measured in 
situ. This data was later used as a quantitative data base against 
which the digitally processed ERTS-1 photographic data were examined. 
Digital Image Analysis 

In order to quantitatively examine image density characteristics, 
digital image processing procedures were developed using TV digitization 
of the Images in registration, storage in a PDP-M computer disc, image 
processing, statistical analysis and display programs. 

Using this digital image processing system, both a spectra! and 
temporal analysis of the reflectance characteristics of the Willcox 
Playa were performed. An analysis of the cloud shadow character! st ics 
was carried out using imagery of Phoenix, Arizona, Imagery of the 
Prescott National Forest was also digitally analyzed. The relative 
effect of the atmosphere i n each of the four spectral bands was mea- 
sured using imagery of the Willcox Playa. The statistics of the rela- 
tive screen brightness range in each of the four ERTS-I bands were 
established for desert areas and a procedure was developed to correct 
for atmospheric effects in each of the four bands. 

2 
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Geologic Interpretation 

A geologic interpretation was performed using ERTS-1 black-and- 
white positives which had been reprocessed to provide optimum additive 
color Image characteristics for exposed soils and rocks found in the 
southern Arizona desert. This reprocessing procedure has been precisely 
set forth herein, A considerable amount of image detail identified by 
color signature differences was rendered visible on the additive Color 
Viewer screen. Many features appearing in the imagery did not appear 
on the most recent geologic maps. 


Image Characteristics Related to Reflectance Spectra 




Section 2 


Techniques and Epuipment 
for ERTS-I Mu I ti spectral Photographic Images 

Mu I ti spectra I Scanner (MSS) imagery exposed on different dates was 
received from the Goddard Space Flight Center in both positive and nega- 
tive form. The spectral bands included the 500-600 nm, 600-700 nm, 700- 
800 nm, and 800-1100 nm spectral bands. Candidate frames from these ERTS 
photographic products were selected, taking into consideration the extent 
of cloud coverage. The selected frames were employed for densitometr ic 
additive color and digital image analysis and correlation with ground 
reflectance spectra. 

Analysis of Photographic Characteristics of the ERTS-I MSS Data 

A photographic analysis of the ERTS-1 MSS photographic data acquir- 
ed over the Willcox test site was performed. The purpose was to establish 
the best procedures for printing positives of semi-arid desert like ter- 
rain such as exists around Willcox, Arizona, 

Release positives of MSS Bands 4, 6 and 7 acquired on 27 September 
1972 were the first Images analyzed. Graphs relating step wedge grey scale 
numbers on the data block to density were generated from the data block. 
These curves are shown in Figure 1. 

The graphs in Figure 1 show that the reproducTicn techniques of 
printing and processing the release positives yielded photographs which 
were well matched for the three bands. However, the relatively high 
minimum density (approximately 0.4) did not allow sufficient image bright- 
ness when the positives were used in an additive color viev/er. The small 
change in density at the higher grey scale values (f.e, for increasing 
relative log exposure) reveals that the release positives were reproduced 

to give greatest density range for the low brightness deta i I In the scene 
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Figure I. Character I stlc curves of ERTS-I release positives of MSS 

imagery of the Wit Icox Playa acquired on 27 September 1972, 

at the expense of the higher brightness detail, 

A densitometr ic evaluation of the reproduction of the release positives 
oi MSS Bands 4, 5, 6 and 7 was performed. Contact prints from the 70mm 
multi spectra I scanner release negatives are shown in Figure 2, Great 
care was taken to match the densities of the step wedges on the prints 
to match those of the release positive transparencies. 

In Band 4 (500-600nm), it is apparent that most of the terrain is 
In highlight. In terms of the step wedge number versus density curve 
(Figure I), the densities lie along the portion of the curve where the 
density differences are not large. The imagery, therefore, appears 
lacking in contrast, especially in the darker cultivated land region 
te the right of the Willcox Playa. The optimization of density dif- 
ferences in the scene, particularly the highlight regions, can be ac- 
complished by Increasing the exposure of the positive image while main- 
taining the development parameters constant. The density range of a 
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500-600 nm Band 


600-700 nm Band 
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700-800 nm Band 


800- II 00 nm Band 


Figure 2. 


ERTS-I bulk multlspoctral scanner Images taken on 27 £eptember 
1972. The Wlllcox Playa Is In the upper right of each photograph. 
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large portion of the scene lies between the and #8 steps of the grey 
scale while the mountain ridges have a lower density corresponding to 
the #13 step. The density differences between most of the terrain ob- 
jects do not exceed ,12 in Band 4. If exposed on the straight line 
portion of the percent transmission versus density curve, the density 
differences produced could have been as large as 3.0. 

A similar examination of Band 5 (600-700 nm) shov/s considerably 
more image detail than is present on Band 4. This increase in detail 
results from an improved contrast between objects in the scene (for 
example, the cultivated area to the right of the Playa). This increased 
contrast does not result from any difference in processing, rather It 
is caused by the fact that the negative material was given sufficient 
exposure to place the objects of interest more nearly on the straight 
line portion of the percent transmission versus density curve. The 
density range of most of the scene lies between steps #9 and #12. The 
image density range is ,3 in Band 5. 

Both near-infrared release positives, Bands 6 and 7 (700-800 nm 
and 800-1 100 nm) were partly exposed in the reproduction, although the 
contrast is the same as Band 5. The high refiectance of the cultivated 
crops adjacent to the Playa can barely be differentiated from the barren 
clay flat of the dry lake bed. This is because both objects, in spite 
of their I nherent spectra I br ightness differences, l ie a long the toe 
of the step wedge number versus density curve. Densitometric measure- 
ments of the Playa image reveals that it has an average density of .52, 
while the crops have a density of ,5. These values correspond to steps 
#8 and #7 on the data bibck step wedge. If the printing exposure had 
been increased, the density differences between these two objects could 
have been increased from ,02 to ,4 in both Band 6 and 7. 
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The set of NASA released positives shown in Fiqure 2 were placed 
Into the Spectral Data Model 64 viewer tor re< onbi nation in additive 
color. Figure 3 Is a photograph of the viewer screen in which Rand 4 
(500-600 nm) was projected as blue, Band 5 (600-700 tm) as qreen, and 
Band 6 (700-800 nm) as red. This renditi r: i- an imat ion of 

that which would be achieved using color infrared film. Because of 
the incorrect printing exposure of the positives, the imagery lacks 
image detail, contrast, and color. 
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Figure 3. Additive color rendition usinq NASA bulk positives in which 
the (500-600 nm) band is projected as blue, the (600-700 nm) 
band as green, and the (700-800 nm) band as red. 


The second set of ERTS mu It I spectre I scanner release positives which 
were analyzed was acquired on 14 October 1972 over the Wlllcox Plava test 
site. The nadir of these photos is cot-'- ' Pfifab ly eastward of that taken 
on September 27th, The relative g^'ev scale number versus density curves 
shown In Figure 4 were generated us inn the grey scale In the data block 
of each release positive. 
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Figure 4. Density vs. step wedqe number of mu It! spectra I scanner 
bulk positives of area covering the Will cox Playa, Arizona, 
taken on 14 October 1972. 

A comparison of this set of curves with those of Figure I shows that 
there exists a larger scene brightness range which Is Imaged on the 
upper portion of Figure 4. In Other words, the exposure has been In- 
creased. There has also been an Increase In the slope of the step wedge 
number versus density curve for these later release positives which are 
shown as contact prints In Figure 5. Each of the multispectral scanner 
positives are we| I matched In curve slope and minimum density. 

I Comparing Band 4 (500-600 nm) of Figures 5 and 2 shows a marked 

Improvement in the Image detail for the cultivated land areas to the 

' _ I. 

soiitheast of the Wlllcox Playa; this. In spite of the relative Increase 

I 

In I density. Is a much better positive Image for analysis In an ad- 
ditive color V I ewer , 
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700 - 800 nm Band 800 - 1 100 ntn Band 


Figure 5, Contact prints of ERTS-I multlspectral scanner bulk Images 
taken on the Wlllcox Playa, Arizona, on 14 October 1972, 
The Playa Is In the upper left of each photo. 




The release positive for Band 5 (600-700 nm) shows that most of 
the scene lies between steps #6 and #10 on the transmission prey scale 
located at the bottom of the Image. This represents a density difference 
of only .3. The most dense images in the photo correspond to step #13 
which lies just at the start of the straight line portion of the percent 
transmission versus density curve shown In Figure 4. A significant por- 
tion of the curve Is not being utilized at all. The highlight areas 
of the WlMcox Playa are In the fog density level of the film so that 
the entire scene exists along that part of the curve which yleid^ the 
least discrimination in terms of density. The base-plus-fog density 
level for all MSS bands exposed on 14 October is excessively high, 
being .49. 

A simulation of color infrared film, shown by Figure 6, was generated 
by projecting the (500-600 nm) band as blue, the (600-700 nm) band as 
green, and both the (700-800 nm) and (800-1100 nm) infrared bands as 
red. The brightness of each near-infrared channel was reduced to half 
value so that when projected their combined brightness would be the same 
as a single Infrared channel projected at full brightness. Figure 7 
shows a false color presentation generated fay projecting the (500-600 nm) 
green band as red, the (600-700 nm) red band as blue, and the (700-800 nm) 
Infrared band as green. 

MSS release positives obtained on 2 November 1972 were typical of 
the photographic characteristics of all subsequent release positives. 

The nominal nadir of the imagery is directly north of the 27 September 
Imagery, Gontact prints of the four release positives are shown In 
Figure 8 and the grey scale number versus density curves associated 
withi them are given in Figure 9. The exposure of the positive images, 
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Figure 6. Composite additive color rendition of 14 October 1972 ERTS-I 
imagery in which the 500-600 nm band is imaged as blue, the 600-700 nm 
band as green, and both the 700-800 nm and 800-1100 nm infrared bands 
as red. 
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Figure 7, Composite additive color image similar to that above in which 
the 500-600 nm band Is imaged as red, the 600-700 nm band is imaged as 
blue, and only the 700-800 nm Infrared band is Imaged as green. 
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Figure 9, Density versus step wedge number curves of ERTS-t bulk 70mm 
positive transparencies of Will cox Playa, Arizona, taken on 
2 November 1972. 

as well as the processing conditions, are almost identical with those 
of 14 October, except that the minimum density is .1 lower. 

The overall density of Band 4 (500-600 nm) I ies between the ^8 and 
0t2 steps of the grey scale which represents a density difference of 
.52. The center of the Playa has an average density of 1.05 and shows 

t 

considerably more contrast within it than either of the previous sets 
of scanner data. Exposure of Band'4 is excessive and no image detail 
lies along the straight line of the step wedge number versus density 
curve. The overall density range of the image detail in Band 6 lies 
between steps #7 and #13 on the positive grey scale. Band 5 appears 
to be much "sharper” than Band 4 in terms of Inherent contrast, although 
the overall density ranges do not differ by more than .28. This Is 
due to a contrast which exists as a result of the spectral reflectance 


of the terrain 


A color Infrared simulation of the 2 November 1972 Imagery has 
been generated by projecting Band 4 (500-600 nm). Band 5 (600-7(X) nm), 
and Band 6 (700-800 nm) as blue, green, and red respectively In a Spec- 
tral Data Model 64 additive color viewer. Figure 10 Is a color repro- 
duction of this rendition created on the viewer screen. The multlspec- 
tral screen brightness factor Is low due to the high base-pl us-foq 
factor. A change In the reflectance of the Playa can be seen by com- 
paring this rendition to Figures 3 and 7. 

Photographic reprocessing techniques were subseguently developed 
to enhance geological Information. These are presented In Section 5. 
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Figure 10. Additive color rendition of 2 November 1972 ERTS-I Imagery 
In which the (500-600 nm) band is nrojected as blue, the 
(600-700 nm) band as green, and the (700-800 nm) band as red. 
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Figure II. Spectral Data Model 64 viewer used for additive color 
analysis of ERTS-I mu 1 1 1 spectra I Images, 


ERTS photographic data products were analyzed using additive color 
viewing analysis technigues. The procudures used, with appropriate 
Instrumentation and significant theories, are described In the subse- 
quent paragraphs. 

Spectral Data Model 64 viewer, shown In Figure II, was an essential 
Instrument for use In the Interpretation and analysis of ERTS MSS Imagery. 
The device produces a single color presentation for the spatially calib- 
rated photos by projecting the Imaqe of one photo on top of the other 
using different color light sources. This technique of analysis permits 
a scientist to select a set of bands and to Interpret within results 
from a single color presentation. In addition, a mul tispectral viewer 
provides the scientist with the capability of altering the color of the 
presentation In order to enhance the particular relationships he may be 
seeking. Fundamentally, the mul tispectral technique allows one to 
create a color presentation specifically for the purpose of his dis- 
cipline and interests. 
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In order to quantitatively evaluate mu I tispectral color Images, 
color must be colorlmetrlcal ly defined as that conscious sensation 
which Is exhibited when light of a specific spectral energy distribution 
enters the eye. It has been experimentally shown that differences in 
this energy distribution cause variations in the observed response of 
the eye and may be described in terms of three distinct psychophysical 
variables. The first is hue which is basically that quality of color 
which leads to the definition of an object as being red, green, yellow, 
etc. Saturation, the second quality of color, is described as the 
amount of white In a given hue. It may be also considered as the con- 
centration of the color. For instance, it is the difference between 
red and pink. As the amount of saturation in a color decreases, it 
approaches pure white. Brightness, which is the third variable of 
color, is described as the amount of visible energy contained in a 
certain hue which is saturated to a specific value. 

DIGITAL IMAGE PRXESSING EQUIPMENT AND TECHNIQUES; 

A digital analysis of the ERTS-I Photograph ic imagery was per- 
formed, The purpose was to make a statistical analysis of ERTS-I pho- 
tographic data in order to determine the spectral and temporal variations 
of the image characteristics, as well as to devise procedures for clas- 
sification of terrain features. 

The equipment used for this purpose was a Spatial Data "System 800" 
This device, shown in Figure 12, and is described in some detail in 

Appendix A, 


Figure 12. The ”System 800” digital process 'ng equipment for 
mu 1 1 1 spectra I photographs. 

This system for digitizing and analyzing mu 1 1 1 spectra I photo- 
graphs operates as follows: 

— Four mu 1 tl spectra I black-and-white transparencies are placed on 

a light box. These photographs are digitized in registration one 
at a time by a TV camera. The transmitted brlnhtness (e.g. the 
density) of each mu I tl spectra I photograph is digitized Into 8 
bit binary numbers (256 density levels) in a 512 x 576 point matrix. 

— The 512 X 576 x 8 bit digitized Image from each mu I ti spectra I 
photograph is stored on the disc of a PDP-ll/E-04 computer. 

~ A series of programs are used to process the images. For instance 
addition, subtraction, multiplication and division of any combination 
of the four multispectral photographs is accomplished by the 
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computer and a series of programs, A description of these programs 
is included in subsequent paragraphs, 

— A digital disc is used to store the processed image and to display 
it on a color TV monitor. The display Is in eight colors, A 
Decwriter also produces hard copy of statistical results and of 
operating program sequence. 

A description of the "System 800" programs is as foi lows: 

MARKER; Reads, saves, and generates joystick- indicated fiducial marks 
for multiple picture registration, 

DIGIT; Digitizes a video picture into a named computer disc 
picture file. 

HISTOR; Generates a line histogram on a black-and-white monitor 

displaying the grey scale distribution of a given mu 1 1 i spectra 1 
picture and accepts joystick inputs to partition the grey 
scale for DISPLAY. 

DISPLAY: Copies a picture from computer disc to the data disc 

using slice values from HISTOR or a uniform default partition. 
AREA; Outputs eight numbers corresponding to the percentages of 
the full picture of each color on the color TV display. 

ADDER: Sets Cpicture 23 = C^picture I) + (picture 2)3/2, where 

pictures I and 2 are named computer disc picture files. 

NEGATE: Produces the compliment of a given picture. 

CONTEX: Performs a contrast expansion of a given picture within 

var iab le I imits. 

SHIFT: Replaces a picture by itself shifted any number of pixels 

up, down left or right. 
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GRID; Generates a variable spaced white-on-black grid on 
the color TV display, 

D I AGON; Performs various exercising, diagnostic and display 
functions under keyboard control , 

Eight bit digitization was selected in order to divide the density 
range of the photographs into 256 levels. The scanner used for digitization 
used a special vidicon-type television camera tube and special circuits to 
accurately sense the light emitted by the image to be digitized. A light 
image of the object is produced on the photosensitive face of the camera tube 
by an optical system containing an ordinary camera lens. Brightness at any 
point on the image is a function of an X and Y coordinate position on the 
tube face. Image brightness is then defined as a third dimension, Z. 

The brightness, or Z-value, is converted into an electrical voltage 
(video signal) by repeatedly scanning the image with an exploring spot 
formed by the electron beam of the camera tube. This spot generates an 
electrical video signal which Indicates the brightness at its instantaneous 
position. The video signal is transmitted through a video amplifier to the 
TV display, where it controls the brightness of a reproducing spot formed 
by the CRT electron beam. The reproducing spot moves over the viewing 
screen in a path similar to, and synchronous with, the exploring spot. The 
reproducing spot reconstructs the brightness distribution in the image area, 
both in magnitude (Z) and position (X-Y) and thus generates a television 
picture of the object. 

Therscanner path covering the image area consists of a series of 
straight parallel 1 1 nes. The spot moves from left to right across each 
line at a constant speed and then returns to traverse the next line. The 


ORIGINAL PAGE B 

OF POOR QUAIJT^ 


20 






spot returns back to the top upon completion of the bottom line. During 


the horizontal and vertical re'^race periods the exploring and reproducing 


spots are disabled, or blanked, so that the image itself is scanned only 


during the forward, left to right motion of the spots. 


In order to reduce flicker in the reproduced picture, the scanning 


lines are Interlaced. That is, instead of moving across the horizontal 


lines in the sequence 0, I, 2, 3, 4, the spot first covers the even- 


numbered lines (0, 2, 4, ..., 478) and then the odd-numbered lines (I, 3, 


5, ..., 479). The scanning beam completes two interlacing scanning 


patterns (fields) to cover the entire image or frame. Each field takes 


1/60 second to complete; the entire image is scanned 30 times a second. 


X ADDRESS 


PIXEL 0 1 2 3 4 


^ V \ 


HOniZONTAL SCAN 
61.2 aS 


V HORIZONTAL RETRAt.'E 

12 3)jS 

VERTICAL RETRACE 

23 LINE PERIODS y 

22 LINE PERlOiiS ' 


Television Scan Pattern 


A sync generator provides the signals for synchronizing the explor- 


ing spot on th#|;'scanner and the reproducing spot on the display. Digital 


counting circuits divide each line into 512 parts while counting the iines 


through each frame. The instantaneous digital values represent X and Y 


positions of the scanning spot. The X-Y addresses select 245,760 points 
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Section 3 


In Situ Spectral Reflectance Measurements 
and Ground Truth Data 

Spectral Reflectance measurements were made In situ at the 
Wlllcox P I aya, Arizona and at Copper Basin and Copper Creek south 
of Prescott, Arizona. The equipment, techniques used and the sig- 
nificant data obtained are described herein. Also described is the 
soil moisture data obtained at Willcox Pfaya and the location of 
homogeneous kinds of trees In Copper Basin. 

Pre ERTS-I Spectral Reflectance measurements of the Willcox Playa 

Willcox Playa, located in the southeastern part of Arizona, is a 

large dry lake which is clearly resolved on ERTS-I Imagery. It is 

bounded on the southeast and northeast by cultivated fields and coarser 

al luvial slopes to the west. The Playa surface is flat and periodic 

fiooding serves as a leveling agent. Variations in sediment size and 

composition, as well as the availability of surface and ground water, 

produce a variety of landforms. The central region of +he Playa is 

hard, compact crust which has no microrelief, except that created by : 

< 

mud shrinkage. The dry surface consists of mud polygons which are a 

' I 

I ight grey buff color. Soft, dry, and porous surfaces are often salt- f 

stained and lead to the sticky-wet areas found where the water table < 

Is near the surface. The soft puffy soils are usually a light brown 

color depending on the salt and moisture content. The Playa is dynamic 

and it changes in response to the surrounding terrain. The surface 

drainage on and adjacent to the dry lake changes as a function of the 

rainfali and a technique is necessary to monitor the dynamics of the 
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area. Spectral reflectance measurements of the dry Playa surface 
Indicate a surprising degree of uniformity. 

Prior to the launch of ERTS-I measurements were made of the in- 
cident solar irradlance spectra at the Wiilcox Playa, Arizona, Simul- 
taneous spectrorad iometr ic radiance measurements were obtained of the 
Playa and large plowed and grass fields adjacent thereto. These data 
were computer processed to provide relative percent directional reflec- 
tance data of homogeneous terrain objects large enough to be well re- 
solved on the space imagery. 

The quantitative spectral analysis of soils, vegetation, and rocks 
requires that simultaneous and accurate measurements of incident and 
reflected radiation be made. Since the spectral energy reflected by an 
object varies with that which is incident upon ft, the spectrorad iometr I c 
measurement must be made at the same instant of time. The arrangement 
of the Instrumentation used to obtain simultaneous measurements of in- 
cident and reflected radiation is shown in Figure 13, Two instruments 
were required to obtain reflected radiance measurements due to photo- 
multiplier sensitivity limitations. 

A spectrorad iometer using a wedge interference filter system enabling 
the spectrum from 380-1250 nm to be continuously scanned was used to mea- 
sure the spectra I distribution of incident solar irradiance. The instru- 
ment was equipped with a diffusing screen so that its directional response 
was proportional to Lambert's cosine law. This technique of measuring 
Incident irradlance was used since solar energy falls upon the earth's 
surface and Is reflected into an entire hemisphere regardless of Its 
original direction of propagation. True cosine response also eliminates 
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Figure 13. Arrangement of Instruments for measuring the spectral dis- 
tribution of Incident and reflected solar radiation in the 380-1250 nm 
region of the spectrum. 


the need for precise aiming to the instrument. 

Whereas a radiometer measures in units of energy rate intensity 
such as microwatts per centimeter square, a spectrorad iometer measures 
In units of energy rate Intensity per bandwidth, such as microwatts per 
centimeter square per nanometer. This latter system of units is most 
meaningful for measurements of radiation since a graph of the spectral 
distribution of radiant intensity yersus wavelength can be obtained. 

The area under such a curve can be made numerically and dimensionally 
equal to energy available In the spectral bands of the ERTS-I multi- 
spectral scanner. The spectrorad Iometer used was capable of measuring 


from 0.01 to 1000 microwatts per centimeter square per nanometer. These 
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values correspond roughly to Fllumination levels of 0.03 to 30,000 
foot candles. 

The calibration of the spectrorad iometer was verified before and 
after the experiment, A spectral standard lamp, serially numbered and 
calibrated against the National Bureau of Standards reference, was used. 
The lamp was of the ribbon filament tungsten type and had a nominal 
accuracy of plus or minus five percent relative to the NBS standard. 

The half bandwidth of the spectrorad iometer used to measure the 
Incident Illumination Is approximately 15 nanometers in the 380-750 nm 
range and 30 nanometers in the 750-1250 nm range. Stray light response 
to unwanted wavelengths of 15 nm bandwidth and far from the v/avelength 
of Interest Is usually in the order of 0,0! percent. 

The periodic calibration of the spectrorad iometer used for measur- 
ing Incident sunlight allowed an accuracy of plus or minus seven percent 
in the long wavelengths and plus or minus ten percent in the short wave- 
lengths of the spectrum. Most of this error, of course, comes from 
uncertainty In the secondary standard used. The relative accuracy of 
all points with respect to each other throughout the wavelength range 
measured Is approximately plus or minus three percent. 

The reflectance radiometer system allowed measurements of the average 
power of the solar radiation reflected by soils and vegetation in the 
Vi/nicox Playa area. By means ot a grating monochromator, these readings 
can be made at selected wavelengths over a bandwidth determined by the 
grating and si its* rThIjs instrument basical ty consists of an optical 
system which limits the entrance of energy to a twelve degree field, 
a monochromator grating to spectrally Isolate the visible energy to a 
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five nanometer ha If band pass, and the infrared ??rierqy tc a ten nanometer 
halfband, as well as detector heads to sense the magnitude of the inci- 
dent energy. 

The reflectance spectrorad iometer is deslqiied so that light re- 
flected from an object passes through a diffuser system and is directed 
by the collective lens into the monochromator housing via the entrance 
slit. Calibration accuracy is obtained when the light incident on the 
diffuser Is imaged on the entrance slit, completely filling the slit 
area with light. A collective lens in the beam input optics in front 
of the monochromator entrance slit, collects the incident ! iqht which 
is properly matched with the diffuser to create a uniform illuralnatlnq 
bundle inside the monochromator housing. This bundle of 1 iqht is then 
incident on a plane diffraction grating where It. is angularly disper- 
sed according to wavelength. Each wavelength present in the source 
bundle reflects off the diffraction greting at a different angle. The 
grating can be rotated to direct any selected v/aveiength bundle onto 
the center of a concave mirror. The mirror collects the light and, 
with the help of a quartz corrector lens, forms an image of the entrance 
slit on the exit si it. 

The visible range grating Is a 1350 qrcxjv© per miUimeter grating 
covering from 350-800 nanometers In the first order and is blazed at 
500 nanometers. The reciprocal linear dispersion is 6,4 nanometers per 
millimeter. The combination of grating and slits determines the dis- 
persion of the systefn. The effective widths of the visible grating was 
20 nanometers as determined by the entrance and exit slit widths are 
2,68 and 1 , 5 m 1 1 I Imeters respectively. These particular grating widths 
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were selected in order to obtain sufficient photomultiplier response 
when measuring targets of iow reflected brightness. 

The visible spectrorad iometer utilized a photo diode detector to 
produce a signal proportional to the intensity of the light which strikes 
It. However, in order to obtain sufficient sensitivity in the infrared, 
a cooled photomultiplier detector must be used to produce an output 

signal which Is of sufficient amplitude to allow significant results 

/ 

to be obta i ned . 

Readout was accomplished using a self-ranging picoammeter. This 
device has the advantage of displaying digital ly the detector current 
output of the spectrorad iometers. Thus calibration and dark current 
values can be readily obtained and monitored. In addition, no scale 
switching Is required and readings can be obtained rapidly* 

In order to obtain accurate data, the spectrorad iometers were 
calibrated before and after the experiment. The calibration of the 
Instrument can shift due to such factors as the collection of dust on 
the optical surfaces or a variety of other random factors. The elect- 
ronic circuitry of the Instrument is very stable and maintains uniform 
response but, nevertheless, the accuracy of the instrument was verified. 
Data Analysis Techniques: 

If a sufficient number of sets of independent measurements of 
Incident and reflected radiation (at 27 points in the spectrum from' 
350-1100 nm) are made and the readings averaged at each wavelength, it 
Is possible to make the Central Limit Theorem apply to the distribution 

of average reflectance values. Thus, 

2-a - C ^ t t- 12£ 
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where; = +he k th measurement of ref leered jadialion at v/ave length 


i, object j, 

y, .. = the k th measuremert of Incident solar radiation at 
I j K 

wavelength ip object ,K 

z. . = the average percent directlcnal reflectance at wavelength 

’ J 

Ip object J* 

n = number of readinus. 

By the Centra! Limit Theorem^ the d i str Ibut jon of z.- wifi be approx- 
imately normal regardless of the distribution of the Individual values of 
Xjji^ and yjjk* This, of course, assumes that measuremenf errors are inde- 
pendent and randomly distributed. 

At each wavelength, for each object, the density spectra is also 
computed : 




The power which would be reflected If the object l;ad Incident upon 
it radiation emitted by a standard sun is also caJculatad as follows: 
Sj. = Zj. X (standard value for each i) 

Both djj and s.^ have been developed to estimate the density of 
the image on multispectra I photograohy. These estimates are currently 
in the process of ref i nement and are auxM iary to the spectral analysis 
which is contained herein. 

A flow diagram of this computer program is shown in Figure 14, 
Note: Computed f or ( 1 ) each wavelength 

• ' ' (2) each object 
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Figure |,4. Flow diagram of percent directional reflectance readings. 

Computer Output: " " 

The computer output format of the computations performed on the 

incidence and reflectance spectra is shown in Figure Associated- 

with each wavelength (first column) is the instrument reading (second 

column) and corrected value of incident sual ight (third column) in 

microwatts per square centimeter. The ref lecta cw measurements and 

* 

associated corrected readings (which are dotermlnod by Hie si'andard 
lamp calibration) are shown in the fourth and fifth columns. 

The percent directional ref tecta nee value shown in the sixth 
column Is the ratio of the' value In the fifth column divided by that 

in the third column. The last column Is the value of a standard 
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DATE: 3/6/69 
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2.35919 
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81.75 
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750 * 

93 

88.9735 
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Figure 15. Computer output format of incident and reflectance spectra. 
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illuminant times the percent reflectance. Thus, the actual power ref- 
lected by an object can be computed by knowing the Incident radiation 
and the data in the sixth column of this output. 

Figure 16 shows the associated computer plot of percent dirsc- 
tfonal reflectance as a function of wavelength. This is an example of 
the percent directional reflectance of grass obtained In situ. The 
chlorophyll reflection band at 550 nm and chlorophyll absorption band 
at 675 nm are clearly evident as is the mesophyil reflection above 
725 nm. 

Incident Solar Radiation: 

The intensity and spectral distribution of solar radiation falling 
upon the earth's surface varies with solar angle and atmospheric condi- 
tions. Measurement of the absolute amount of solar energy which strikes 
the terrain Is important since all ground objects reflect different 
amounts of radiation in each spectral band which is in direct propor- 
tion to that which is incident. 

Theoretical analysis of solar Illumination based upon air mass cal- 
culations and Raleigh scattering are of little practical value in pre- 
dicting the spectral distribution of solar energy wh ich actua ( ! y reaches 
the ground. This is primarily due to the existence of unknown amounts 
of Mie scattering and absorption In the atmosphere due to particles 
which are large compared to the wavelength of the radiation. 


Figure 17 demonstrates this condition. Spectroradiometric readings 
of solar energy (both suni Ight and diffuse skylight) using a jambertian 
detector were measured during a four day Interval at WlilcoK, Arizona, 
Spectral intensity in microwatts per centimeter squared per nanometer 
from 380-1200 nanometers is shown for four of these measurements between 
1104 and 1515 MST on March 6. 





WmSx'pfaya^tesrsite^^^ °l incident at the 

riaya rest site at different times during the day. 


The first feature cleariy evident from this data is the iarge varl 
intensity with wav/eiength as a function of time (solar angle). 
The characteristic absorption band beiow 380 nanometers at 725 and 950 
nanometers is also clearly evident. 

The percent directional reflectance of the till Icox Playa, Arizona, 

was obtained for two adjacent areas at the location shown in Figure 18 . 

Simultaneous apectrcradiometric measurements of Incident global solar 

radiatloni were made along with measurements of the radiation reflected 
by the Playa. 

The reflectance spectroradlometers were oriented so as to prevent 
any specular reflection from the Playa surface from ehterlng the input 



SPECTRA OF PLAYA 


TARGET ARRAY 


n^nr COCHISE, ARIZ. 

OK N»00— W1094S/1S 


CONTOUR INTERVAL 2» FEET 
DATUM •• MtAM UA kCVlL 


Figure J8. Location at which Incident and reflectance spectra were measured. 








Figure 20. Percent directional 
reflectance of bare soil, crass, 
and the Willcox Playa. Note 
the comparative spectral uni- 
formity of the Playa compared 
to higher infrared reflectance 
of grass and bare soil. 


-n N, 


I • 



\ 


optics of the Instrument, The partly cloudy atmospheric conditions 
that existed requlr^ that extreme care and patience be used In making 
these measurements. The Incident and reflected spectrorad lometrlc data 
were corrected and the percent directional reflectance calculated using 
the computer program discussed previously. 

The percent directional reflectance data presented In Figure 19 
Indicated a surprising degree of uniformity In the reflectivity of the 

i 

I 

Pilaya surface. Although there appears to be slight brightness differ- 
ences between the two areas, the shape of the reflectance curves Is 
quite siml lar. 

A comparison of the spectra of plowed fields and grass with the 
Playa Is shown In Figure 20, The area from which these readings were 
obtained was adjacent to the target array at the locations shown In 
Figure 18, The characterlst Ic chlorophyll absorption band (675 nm), 

chlorophyll reflection band (550 nm), and mesophyll reflection band 

! 

(bbove 700 nm) are shown In the spectra of the grass. 

Additional Reflectance Spectra of Willcox Playa 

The reflectance spectra of twelve locations In the Willcox Playa 

I 

was obtained In the Fall of |973, The precise locations recorded on 
two separate occasions are shown In Figure 22, The appearance of the 
Playa from the roadway Is shown In Figure 23; Figure 21 Isa photograph 
of the Playa also looking eastward. 

The incident solar radiation falling on the Playa measured In 
micro watts per square centimeter per nanometer at f0»22 POT on 4 Septem- 
ber 1973 is shown In Figure 27, 
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Figure 21. Willcox Playa, looking east. 

The percent directional reflectance of the twelve locations on 
the Playa is shown in Figures 28 and 29. Considerable small local 
differences were found to exist about site I due to soil moisture 
variations and hence three sets of reflectance spectra were taken at 
this location. These variations in the apparent condition of the 
Playa are shown In the photographs In Figures 24 through 26. 

Table I, lists the percent soil moisture at each location. The 
soil samples were taken at the same time as was the reflectance spectra. 
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I Figure 22, Locations In Wjllcox Playa, Arizona, where Spectra was taken during Fall 1973. 



Figure 23. Wi I Icox Playa, viewed from road looking east. 



8.7710 

9.5476 

10.6016 

8.6695 

8.2077 

7.0687 

7.7994 

8.9289 

8.1270 


I I .4335 


8.8666 


9.9244 

8.7659 

9.0585 


5.851 
5.486 

67.672 

4.191 

3.708 

3.851 
4.417 
4.288 
3.559 


3.428 


3.068 

4.590 

3.392 

3.663 


Table |. Soil Moisture at different locations In Wlllcox Playa, Arizona, 
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Willcox Playa location I A, soil moisture 5.9 percent 


Figure 24 




Figure 26. Willcox Playa location 
1C, soil moisture 67.7 percent 


Figure 25. Willcox Playa location 
IB, soil moisture 5.5 percent 
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Figure 27, Incident Sun Spectra taken in Wlllcox Playaon 4 September 1973 
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Figure 28. Percent directional reflectance at different locations In WIIIcox Playa on 4 Sept. 1973 
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at different locations in 
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The percent soli molslure data obtained were plotted against the 
percent reflectance at 650 nm and 950 nm shown below. 


Percent directional reflectance 

at 650 nm vs. percent soil moisture 

at various locations in Wlllcox Playa. ® 


•/o Soil Moist uec 


Figure 30, Percent reflectance at 650 nm vs, percent soil moisture. 
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Percent directional reflectance 
at 950 nm. vs. percent soli moisture 
at various locations In Wlllcox Playa, 


^ % SOIL Moisture 

Figure 31, Percent reflectance at 950 nm vs, percent soil moisture. 
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Figure 32. Incident Sun Spectra taken in Willcox Playa on 16 November 1973 
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As can be seen from the two graphs In Figures 30, and 31, soli 
moisture and percent reflectance at these two wavelengths are un- 
correlated on the Wlllcox Playa. Since water absorbs Infrared radia- 
tion, the effect of trace amounts of water moisture on reflectance was 
to be analyzed spectrally. 

Figure 32 Is a graph of the Incident solar Irradlance at the 
WlMcox Playa measured In micro watts per square centimeter per nano- 
meter, The data was taken on 16 November 1973, at 1020 POT, Nine of 
the twelve locations shown In Figure 22 were recorded and the percent 
directional reflectance measured. These data are shown In Figures 33 
and 34, 

Reflectance Spectra of Predominant Vegetation In the Prescott 

National Forest 

In situ spectral reflectance measurements were taken of three pre- 
valent vegetation types In the Prescott National Forest, These species 
are Pinion Pine, Juniper, and Manzanita, Typical reflectance spectra 
are shown In Figure 35. 

Reflectance spectra In Figure 35 were obtained between 23 October 
and 28 October 1972, 

A computerized telespectrorad lometer was utilized to collect the 
In situ reflectance spectra. This truck-mounted Instrument Is positioned 
adjacent to the target to be measured (e.g., a Pinion Pine) and the de- 
tector head extended In the air to obtain the similar downward perspec- 
tive as the ERTS mul t Ispectra I scanner. The Instrument measures and 
stores In the computer memory the reflectance of a standard white plate 
which Is placed In the field of view of the spectrorad lometer. The 
reference plate Is then removed and the spectra of the target Is obtained. 
The computer enters all the necessary Instrument corrections, divides 
the standard white plate reference readings Into the target readings at 
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Figure 36. Truck-mounted telespectrorad lometer. 

each discrete wavelength measured and prints out the result on a tele- 
typewriter. Figure 36 shows this truck-mounted Instrument In operation 
at the Prescott test site. 

The detector head mounted atop the extended ladder Is comprised 
of an imagery optical system, a circular Interference filter monochro- 
mator, visible and infrared photomultiplier detectors, and an electronic 
amplifier. The circular filter is Indexed Into position In the optical 
svstem In front of the detector which selects the particular wavelength 
to be measured. Usually 80 discrete wavelengths are measured between 
400 nm and I 100 nm. 

The ratio of target reflectivity to standard reflectance reference 
is automatically calculated by the Instrument and displayed under computer 




control. The computer is programmed with a list of spectral ranges 
(visible and Infrared), data points (discrete wavelengths of interest), 
and the number of scans per data point. The computer Incorporated in 
the spectroradlometer Is a 16-blt word length general purpose computer 
having a 4096 work (8K byte) memory. Four hardware accumulators are 
used in addition to a PC register. The basic control program resides 
in| 3000 core locations and inciudes a floating point interpretive package 
The 1000 remaining locations are used for temporary data storage during 
target data acquisition. The basic control program takes radiometric 
measurements at pre-determ i ned wavelengths, stored in a look-up table. 

For reflectance measurements a white standard reflector is compared to 
the target and the relative reflectance of the target is computed as a 
ratio of the target to the standard. The computer then uses this res- 
ponsivity data to calculate the radiance of a target. 

The telespectroradiometer used to collect the data contains two 
spectral channels - one covering from 400 to 700 nm (visible), the 
other from 675 to 1350 nm (infrared). The system is comprised of an 
optical head which contains the telescope optics, a sensing head con- 
taining the photomultiplier detectors, and input-output control console 
which contains the computer, motor drives, ana log-to-digita I converter, 
and teletypewriter. 

The system is designed to have flexibility in the spectral region 
of operation and in its mode of operation. The primary filter assembly 
of| the system utilizes a circular variable Interference filter which 
allows continuous variance pf the spectral region transmitted to the 
detector as a function of the angular position of the filter. The 
spectral response typically covers a two-octave range of wavelengths 
with an average spectra I reso I ut ion of 2 1 12 %, The system operates i n 


an Incremental ly stepped mode which provides the ability to take precise 
measurements at each wavelength. Time Is allowed for gain ranging, amp- 
lifier Sfett ling, a nd multiple read I ngs for computer averag I ng , 

The automatic telespectroradlometer used to collect the data con- 
sists of four basic units: the optical head, the sensing head, the con- 

trol cabinet, and the output unit. 

The telespectroradlometer contains a 300mm focal length refractive 
optical system with a 50mm diameter entrance aperture. The mechanical 
configuration of the optical head of the system Is shown In Figure 37, 
The Incoming radiation Is focused by the objective lens on an Inter- 
mediate focal plane where an adjustable field stop Is located to define 
the system’s field of view. 



Figure 37, Telespectroradlometer optical schematic. 
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The radiation passing through the field stop Is refocused by the 
transfer optics onto the front surface of the circular variable filter. 

The radiation traversing the filter Is refocused by a relay lens onto 
the detector. 

The circular variable filter is motor-driven under computer control. 
The direct coupled angular encoder provides a feedback signal to the 
computer of the instantaneous position of the circular variable filter. 

The computer program specifies the angular position at which measurements 
will be made and monitors the opt i car encoder output until the correct 
angular orientation is achieved. 

The electrical signals generated by the detector(s) are amplified 
by a picoammeter amplifier with remote gain-ranging controlled by the 
computer program. The high-level, low-impedance analog output of the 
amplifier may be transmitted up to a distance of 250 feet to the computer, 
allowing completely remote operation. Also located In the optics package 
are the shutter and viewing mirror solenoids, optical encoder buffer 
amplifiers, and line drivers. 

The picoammeter is a two-stage, solid-state amplifier using compact 
metal-film feedback resistors and solid-state switching for range selec- 
tion, The program-controlled or manually-controlled shutter allows dark 
current measurement and corrections to be made under program control with- 
out operator intervention* A block diagram of the electronic circuitry 
Is given In Figure 38. 

The analog signal from the optical detector Is fed to the analog 
multiplexer of the analog-to-digital Input channel. The multiplexer 
selects the picoammeter output for conversion or other accessory func- 
tions, such as internal self-check functions. The multiplexer output 
drives the sample-and-hold and the eleven-bit analog-to-dlglta I converter. 



Figure 38, Telespectroradlometer electronics block diagram. 

The multiplexer, sample-and-hold and analog-to-dlglta I converter are 
controlled by the MADC Interface wh I ch a J so conta I ns the picoammeter 
gain control registers, shutter and viewing mirror solenoid buffer 
registers. 

The encoder outputs drive the decoder module which. In turn, 

.drives the present position register In the CVFC Interface, The present 
position register supplies the computer program with the current angular 
rotational position of the circular Interference f I Iter, The GVFC In- 
terface contains auxiliary buffer registers for operating the motor 
drives and synchronization control modules. 















SECTION 4 


DIGITAL IMAGE PROCESSING OF ERTS-I PHOTOGRAPHIC DATA 

A digital analysis of ERTS-I Photographic data products was per- 
formed using the computerized "System 800”. This system digitizes multi- 
spectral photographic transparencies on to the disc of a PDP-11 computer; 
processes the data using a series of image analysis programs; displays 
the resultant data on a color TV monitor and prints out the associated 
hard copy on a Decwriter. The description of the "System 800” is con- 
tained in Section 2 and in Appendix A. 

Analysis of ERTS-I imagery were performed using the "System 800” 
of the fol lowing; 

~ Will cox Playa Arizona, Spectral image analysis 

~ Will cox Playa Arizona, Temporal image analysis 

— Phoenix Arizona, Spectral image analysis 

— Prescott Arizona, Spectral image analysis 

— Analysis of atmospheric effects on ERTS-I image characteristics 
Willcox Playa Spectral Image Analysis 

A spectral image analysis was performed of NASA ERTS-I release 
negatives dated 5 June 1973, of the Willcox Playa using the "System 
800”, The purpose was to investigate the spectral uniformity of the 
playa surface on that date. 

Bands 4, 5, 6 and 7 and a reference grey scale were placed on the 
digitizer I ight table at the same time. The digitizer maximum and 
minimum vid^o signals using a linear digitization were adjusted so 
that the playa and surrounding area in Band 5 fall wjthin the 256 level 
digitization capabilities of the instrument. This maximum digitization 


level corresponded to grey step number 4 and the mlnimurif! level to 
grey step number 15, A histogram of the 256 levels of the 250,000 
pixels digitized from band 4 were examined and bands 4, 5 and 6 
ware digitized as foilowes; 


adjusted as follows: 


0-28 

= 

Black 

29-43 

= 

Red 
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0-30 

= 

Black 

31 750 
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Red 
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Green 
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Yellow 

91*110 
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Blue 
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Magenta 
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Cyan 
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White 


Figures 39 and 40 show a contact positive of each input ERTS 
negative and the resulting 8 color display of the digitized Image 
produced by the "System 800" of the Will cox PI aya portion of each 
Input ERTS feh Ip. Since negatives were Inputed into the "Syst^ 800" 
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Band 4 (500-600 nm) 


Band 5 (600-700 nm) 


9, WILLCOX PLAYA "SYSTEM 800" DIGITAL IMAGE ANALYSIS OF 
5 JUNE 1973. Band 4 and Band 5 negatives. 


Original 








Band 6 (700-800 nm) 
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Band 7 (800-1100 nm) 


Figure 40. WfllcoxPlaya "System 800" Digital Image Analysis of 
ERTS-I 5 June 1973, Band 6 and Band 7 negatives. 







brightness is inversely related to the digitization levels and colors 
listed above. 

The Playa is brightest in Band 4 where the 0 to 73 digitization 
levels contain all Image detail with at least 60 percent of the Playa 

falling Into the 59-73 level. The extreme brightness of the extreme 

! 

south west portion of the Play© Is evident. Band 5 shows lower bright- 
ness than does Band 4 with about 50 percent of the Playa image detail 

i 

falling Into the 44-58 brightness category without any of the extreme 
bffightness detail in the south west part of the Playa that was evident 

. i , ' 

iiji Band 4. Band 6 is less bright than Band 5 In all areas of the 

i 

pjaya although the south west area is still brightest. Band 7 contained 
the least brightness of all the four ERTS-I chips although the difference 
In reflectance across the surface of the Playa was similar in all 
four ERTS-I bands. 

Willcox Playa - Detailed Spectral Image Analysis 

Positive transparencies of the four ERTS-I spectral bands taken on 
4 September 1973 were digitized in registration using the "System 8O0^ 
¥(ith MARKER, HISTOGRAM, DIGITIZE and DISPLAY programs. The size of the 
image digitized was made as large as possible within the limits of the 
video camera optics. 

The digitized positives of Willcox Playa for each of the four 

I 

ERTS-I spectral bands are shown In Figure 41. Superimposed on the 
Imagery Is the frequency distribution of the 294,912 digitized pixels 
in each band. 
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Figure 41, Digitized ERTS-I 4 September 1973, Images of Wlllcox 


Playa from positive fllw transparencies showing statistical distribu 


tion of brightness levels (dark at top of Vertical scale In photos) 
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The AREA Program was run on each of the four bands using the 
digitized grey levels shown In the following table. 


D i g i 1 1 zed 
Grey Levels 

4 September, 1973 

ERTS-I Positive 

Band 4 

Band 5 

Band 6 

Band 7 

^ 0-32 

4.32 

3.12 

5.64 

2.38 

33-64 

4.70 

1.98 

0.62 

1 .00 

65-96 

10.02 

4.30 

1.36 

2,24 

97-128 

26.30 

10.04 

5.90 

5.98 

1 23- ( 60 

37.18 

26.30 

27.62 

39.94 

161-192 

16.26 

35.76 

46.38 

45.90 

1 93-224 

1.04 

17.62 

12.00 

2.32 

225-255 

0.08 

0.76 

0.36 

0, 14 

Totals 

99.9 

99.88 

99.88 

99.9 


Table 2. Digital Brightness Classification of ERTS-I 
mujtlspectral images of Willcox Playa taken on 4 September, 1973. 

A histogram has been plotted of the digitized brightness levels 
ip each of the four ERtS-l spectral bands and have been plotted In 
figure 42. Table 3 below shows the expected value and standard 
deviation of the brightness in each of the four bands. 



4 September 1973 ERTS-I Positives 1 


Band 4 

Band 5 

Band 6 

Band 7 

Expected Value, E (x) 

124.9 

155.4 

156. 1 

152.8 

Standard Deviation, { 

41.2 

45.0 

44.4 

34.0 


table 3, Statistics of the 4 September 1973 ERTS-1 positives of the 
Willcox Playa, Arizona. 
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From the above statistics, the histograms shown in Figure 42 

and the frequency distributions superimposed on the ERTS-I digitized 

images In Figure 41^ the following Inferences are made; 

The Playa Is darker In Band 4 than In any of the other 
ERTS-I Bands. 

— The brightness of the Playa Is approximately the same In 
Bands 5, 6 and 7, 

~ The variability of the brightness data is less In Band 7 
then In any of the other ERTS-I bands. 

— The frequency distribution of brightness has two modes 

(peaks) in bands 4 and 5 which are absent from bands 6 and 7. 


Wlllcox Playa - Temporal Image Analysis 

In order to investigate the change In reflectance of the Wlllcox 
Playa as a function of time, an analysis of four ERTS negatives obtained 
on four different dates were made. 

The following four ERTS-I negatives, a 1 1 Band 5 , were used as 
Inputs to the "System 800" and digitized In the linear mode: 

I November, 1972 
18 February, 1973 
18 May, 1973 
5 June, 1973 

All four bands were digitized in registration using DIGIT and 
MARKER PROGRAMS. Following this, the HISTOGRAM Program was run class- 
si fying the 256 digital levels Into exactly the same light equal levels 
on all of the multidata Imagery. The AREA program was then run and 

I 

used to compute the percent of the 294*912 picture elements falling 
Into each of the eight classes. The results are shown In the fable 4 
on page 65, 
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I November i 972 


18 February 1973 


original 





grey levels 
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mmm 


! 0-32 
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1 'i 

Black 

Ii.l8 

14,48 

1 

.00 

3.28 

1 33-64 

1 

, Red 

2.72 

i 

1.98 1 

0 

.58 

31.26 

65-96 

Green 

B.40 

8.70 

6 

,90 

48.46 

1 97-128 

Yel low 

37.90 

19.92 

20.32 

14.06 

129-160 

Blue 

29.24 

35.06 

42.92 

2.52 

161-192 

i 

Magenta 

li.76 

13.26 

18.88 

0.32 

193-224 

Cyan 

5.10 

3.76 

f 

).68 

0 

225-255 

White 

3.62 

2.72 


!.62 

0 


Tota 1 s 

99.92 

99.88 

99.90 

99.90 


Table 4, Digital brightness classification of mu I tldate Imagery of Wlllcox Playa. 


The display program was run and each brightness classified band 
was photographed and Is reproduced In Figure 43. In addition, the 

18 February 1973 and the 18 May 1973 Band 5 negatives were digitally 

. I 

avjeraged, A visual comparison of this result shown In Figure 44 

with the individual brightness classifications for these dates depicted 

i. , , , 

In Figure 43, shows very little variance between these four dates 
for the Playa Itself. 

A histogram was constructed of the digital brightness clas- 
sification of Band 5 on each of the four dates which Is shwn In 
Figure 45. The mean of this data Is tabulated below. 


^nd 5 - Date 

Mean (x) 

1 November 1972 

I26«7 

18 February 1973 

120.7 

18 May 1973 

143.9 

5 June 1973 

74.2 
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18 February 1973 
ERTS-I Band 5 


PAGE IS 
quality 


origwal 
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Average of 18 February and 18 May 1973 
ERTS-I Band 5 


Figure 44. Digitally averaged brightness values of the Wlllcox 
Playa, Arizona on two dates. 
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Figure 46 shows the percent of each of the brightness levels on the 
four multidate negatives in each of the eight digital classes. Ciear~ 
ly the 5 June 1973 imagery is anomalous in each class to the other dates. 

Band 4 of four ERTS-I positives of Willcox Playa obtained on the 
four different dates were digitized in the log mode. These four digit- 
ized Images with superimposed frequency d istributlons of brightness 
levels are shown in Figure 47. The digitized grey levels as computed 

using the AREA Program as shown In Table 5. A Histogram of these bright- 
ness levels is plotted in Figure 48. 
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J i g 1 1 1 zed 
grey levels 

Band 

4 ERTS-I Imagery of Willcox Playa I 

4 September 1973 

18 May 1973 

n 

November 1972 

18 February 1973 

0-32 

3.12 

5.50 


0.96 

20.92 

33-64 

i.98 

5.92 


0,08 

24 L 4 O 

65-96 1 

4.30 

21.84 


0.36 

24.72 

97-128 

10.04 

30.04 

I.IO 

15.26 

129-160 

26.30 

19.70 

1 1 .38 

10.72 

161-192 

35.76 

14.72 

64.70 

3.60 

1 93-224 

17.62 

1 .98 

18.44 

0.28 

225-255 

0.76 

0.22 

2.82 

0 

Tota 1 s 

99.88 

99.92 

99.84 

99.90 


Table 5. Digital brightness classification of multidate Imagery of Wtllcox Playa. 

In order to evaluate the change |n reflectance of the same areas of the 
Playa on the four different dates, three areas of the Willcox Playa were 
measured in each photograph: I) the bright area at the south end of the 

Playa, 2) the west area of the Playa adjacent to the railroad track, 3) lo- 
cations within the body of the Playa chosen randomly. The results are 
presented in the table on Page 72, 
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18 February 1973 Band 4 


I November 1972 Band 4 


IB May 1973 Band 4 4 September 1973 Band 4 

Figure 47. Digitized ERTS-I Band 4 Images from positive film transparencies 
showing variations (.5 to 


6y) reflectance on four different takes 
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I 
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Brightness measurements In the Playa depend on location. The 
higher south end of the Playa Is brighter and shows less deviation 
from the mean brightness than the main body of the Playa. The area 
of the Playa north west of the railroad tracks crossing the Playa, 
where water Is known to be Impounded, shows the greatest variability 
on the four dates. 

From these dates It Is clear that variations in the spectral 
reflectance of ‘the Witicox Playa as recorded by ERTS-I vary greatly 
between spectral bands and In the same spectral band on different 
dates. 


Phoenix, Arizona Spectral Analysis 

A spectral analysis of GPTS-I Imagery of an urban area of Phoenix, 
Arizona taken on 14 April 1973 was digitally analyzed to establish atmos- 
pheric effects on the Image in urban scenes in seml-arld terrain. Four 
carefully prepared positive Images were made from NASA release ERTS-I nega- 

1-lves. These were I inearly digitized in registration using the MARKER and 

i 

DIGIT program; following which HISTOGRAM and AREA programs were run using 
about equgi digitization levels as shown in Table 8 on the next page. 
Digital brightness measurements were made of bright cloud tops 
and dark ground shadows which were cast over urban areas of Phoenix. 


The mean and standard deviation of this data are tabulated in Table 7, 



Cloud Tops 

Shadows 

mean ~ 

• Standard 
deviation (s) 

mean x 

Standard 
deviation (s) 

Band 4 

170 

20 

52 

12 

Band 5 

173 

19 

46 

12 

Band 6 

210 

7 

69 

12 

Band! 7 

223 

22 

63 

12 


Table 7, 




Digital brightness statistics of clouds and shadows appearing 

In Phoenix, Arizona 14 Apr 

,75 ■ 




0-32 
33-64 
65-96 
97-128 
129-160 
I6|)-I92 
I ^3-224 

I 

225-255 


Table 8. 


1 

4 April 1973 

ERTS-I Image 

1 IsSQEJKIHlii 1 

13.56 

I5.!20 

! 1.96 

2.52 

21.70 

|l6j66 

i 4.54 

5,62 

26.42 

23.86 

i 12.62 

15.44 

30.92 

26.76 

24.80 

30.68 

6.64 

14.96 

31.82 

28.76 

0.66 

2.18 

20.60 

13.48 

0.02 

0.30 

3.50 

2.94 

0 

0 

0.08 

0.42 

99.92 

99.92 

99.92 

99.86 


Color 

Code 


Black 

I 

Red 

Green 

Yellow 

Blue 

Magenta 

Cyan 

White 


Totals 


Digital brightness classlf Icatlon of urban section of ERTS-I 
Imagery of Phoenix, Arizona taksn oh 14 April 1973, 


From the above data It Is seen that the difference In the means of 
the brightness of cloud tops and shadows increases as follows; 


Band 4 
Band 5 
Band 6 
Band 7 


mean (>^‘*>»Xt) 
1 18 
127 
141 
160 


Standard ^ 
deviation (s*4- st) 


Table 9, Brightness difference between cloud tops (x, ) and shadows 
(x^) In Phoenix, Arizona 14 April 1973. 

Assuming comparatively uniform spectral reflectance for cloud tops 
and cloud shadows (In an urban area), we see that the contrast In each 
band as a percent of Band 7 Increases with Increasing wavelength as 


fol lows. 


Band 4 _ ^ 

Band 7 ^ ! 00 = 74? 


Band 6 


X 1 00 = 88? 










Prescott, Arizona Spectral Analysis 

Four ERTS-I negatives taken on 29 September 1972 were digitally analyzed 
In order to establish if known communities of Ponderosa Pine could be detected 
in the Copper Basin Area of the Prescott National Forest. The results of the 

HISTOGRAM and AREA programs run with the digitization levels (partitions) 

\ 

are as fol lows: 


Digitized 
grey levels 

0-32 

33-64 

I 

65-96 
I 97-128 
129-160 
161-192 
1 93-224 
225-255 


September 


Sgatives 


Code 

Band 4 

Band 5 

I Band 6 

Band 7 

Black 

4.34 

9.64 

j 9.22 

7.00 

Red 

11.96 

16.80 

24.84 

23.38 

Green 

27.54 

25.84 

33.16 

36.10 

Yellow 

29.34 

19.28 

18,86 

20.06 

Blue 

15.50 

11.82 

9.30 

8.56 

Magenta 

7.18 

6i92 

2.90 

2.52 

Cyan 

3.16 

4.00 

0.78 

0.84 

White 

0,90 

5.60 

0.86 

1.44 

Tota 1 s 

99.92 

99.90 

99.92 

99.90 


Table 10, Digital brightness classification of Prescott National Forest 
on ERTS-I negatives taken 29 September 1972. 


Positive prints of portions of each of the ERTS-I bands are 

I . ■ 

reproduced In Figure 50. Textural patterns in Copper Basin Area 
are evident In this imagery principally in Band 7, The results of 
the "System 800" are shown in Figure 51. These data are analyzed 
with respect to ground reflectance spectra in Section 3 of this report. 
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50. Prescott, Arizona - Positive ERTS-I Imaqes acquired 
on 29 September 1972. 
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Average of Bands 4, 5, 6 and 7 

Prescott, Arizona - digitized brightness presentation for 
photographic negatives acquired by ERTS-I on 29 September 1972 







Atmospheric Effects In Desert Areas 

The relative effect of the atmosphere In each of the four spect- 

» 

ral bands was anallzed using ERTS-I Imagery of the W1 Mcox, Arizona, 
area acquired on 4 September 1973. Positive photographic Images of a 
desert area containing clouds were log digitized. The digitization 

range was adjusted to give it the same range for all four photos in 

1 

order to maintain a constant relative relationship of the eight digit- 
ization levels between the four spectral bands. 

By use of the MARKER, DIGIT, HISTOGRAM, AREA and DISPLAY programs 
the following brightness classifications of each image were made: 


D ijg i t i zed 
grey levels • 

1 4 September (9/3 ERIS-I Positive ! 

BanHT^ 

Band 5 

Ba nd 6 

Band / 

0-32 

1 .30 

1 1.56 

1.02 

1 .24 

33-64 

0.46 

1.80 

0.68 

0.58 

65-96 

1.62 

4.80 

1.34 

1 ,30 

97-128 

9.64 

9,42 

3.74 

3.04 

1 29- 1 60 

62.80 

49.28 

38.46 

24.78 

161-192 

20.92 

30,34 

51.48 

65,38 

193-224 

2.34 

2.14 

2.62 

2.84 

225-255 

0.92 

0.58 

0.54 

0.74 

Totals 

■■ ■9'9".9'D 

99.92 

99.88 



Table U, Digital brightness of ERTS-I images of Southern Arizona 
taken on 4 September 1972, 

j . Ten digital brightness measurements of clouds, cloud shadows 
and desert terrain immediately adjacent to cloud shadows were made 

f 

using HISTOR program. The mean (x) and standard deviation (s) of 

■I ■ ■ 

! , ■ 

these measurements are tabulated on the following page Table 1 2 

The digitization levels were adjusted to make the brightness of cloud 
approximately the same in ell bands. 
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4 Se 

otember 1973 E 

RTS-I Positives | 

Image Area 
measured 

_ Band 4 
“x s 

Band 5 

X s 

Band 6 
X s 

Band / 

X s 

Bright Cloud 

230.7 10.4 

230.3 4,3 

229.0 1.8 

232.9 2.4 

1 

Terrain in shadow 

99.4 7.7 

59.6 12.9 

51.7 8.6 

20.1 10.8 

Terrain not 
In shadow 

145.5 3.6 

141. 1 3.7 

154.1 2.2 

156.4 4.6 


Table 12. The mean and standard deviation of Image brightness measure- 
ments of desert terrain in shadow and illuminated by sunlight. 


The full printout of all measurements, the statistics of which are 
shown above is contained in Appendix F . The statistics in this table 
were computed from a sample of ten (e.g. n=IO) random measurements of 
each of the three areas on each of the four ERTS-I positive images as 
fol lows: 

•vt 

sample mean = x = . 

i y 

sample standard deviation = s = 

A plot of the digital brightness levels of all four ERTS-I bands 
compared to band 7 as a standard Is shown in Figure 52. The increas- 
ing brightness of the terrain in shadow with decreasing brightness and 
decreasing wavelength is evident from this graph. This data quantifies 
the effect of atmosphere with brightness and wavelength. 

j Correction factors were computed to make the digital brightness 

, j • _ 

values in all four bands equal to the same value (e.g. 230). These 
corrected values are shown in Table 13. 

The percent decrease in scene brightness range with decreasing 
wavetength from the last row In Table 13 has been plotted in Figure 53 , 
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Ffgure 52, Digital brightness levels of the four ERTS-I bands compared to 
Band 7 for cloud tops, terrain in sun, and terrain in shadow. 



4 September 1973 ERTS-I pos 

itives 1 

■sEDEKHi 

■SIsfli 

bBand 6 


Correction factor 

.9969657 

.9986973 

1 .0043668 

.9875483 

Bright cloud x corr. factor 

230 

230 

230 

230 

Terrain in shadow x corr, factor 

99.1 

59.5 

51,9 

19.8 

SunI Ight terrain x corr, factor 

145.1 

140.9 

154.8 

154.5 

A; = corr. factor x (sunlight-shadow! 

46.0 

81 .4 

102.9 

134.7 

A as a percent of Band 7 A 

34.1 

60.4 

76.4 

100 



Table 13 , Terrain and shadow measurements normal i zed to equal cloud 
brightness measurements, 
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Figure 53* 


Wavelength (nanometers) 

Scene brightness range of sun I It desert terrain minus 
similar terrain In shadow. The range of Band 7 Is used 
as 100 percent. 


* 

4 September 1973 ERTS-I positives 

dand 4 

Band 5 

■sgnpgf 

ISH 

I (19.8) (corr. factor) t 

•r (terrain In shadow) 

.20 

.33 

.38 

1.00 

2 (154.5) (corr. factor) t 

T (terrain in sunlight) 

1.06 

I.IO 

1 .00 

1 .00 

3 reclprlcal of row 2 above 

.943 

.909 

1.00 

1.00 

row 1 times row 3 ~ corr. facte 
to normalize scene brightness 

r 

.19 

.30 

.38 

1. 00 


Table 14. (]lomputatlonal steps to arrive at correction factors to compen- 
sate for atmospheric scattering of low brightness desert terrain. 



In order to derive correction factors which can be applied to 
each band to compensate for atmospheric scattering as a function of 
wavelength, all data were normalized with respect to Band 7 which 
is used as the standard since It Is known that Rajelgh scattering Is 
least In the Infra-red wavelengths. 

The data from Table 14 have been used to construct a graph 
showing the approximate correction to apply to ERTS-I positive Images 
of desert areas as a function of brightness In each band. The pro- 
cedure Is: 

1 - Adjust digitization so that bright cloud tops have equal 

brightness values in all four bands, 

2 - All image areas in deep shadow should have digitized bright- 

ness values greater than zero, 

3 - Measure adjacent desert area in shadow and in sunlight in all 

ERTS-I bands. The brightest value = 100 and the darkest 

• value = 0 . 

4 - Apply the correction factor shown in Figure 54, 

A summary of the analysis of variability of the digital brightness 
measurements in each ERTS-I spectral band is shown In Figure 55, 

j 

Tliereln It Is seen that the statistical stability of the variability 

i 


of the measurements of digital brightness of sunlit and shadow desert 
terrain does not vary with ERTS-I band. However, the variability in 
cloud brightness measurements decreases with Increasing wavelengths. 
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Section 5 


1 




ERTS Geologic Interpretation using Optimum Image Character Istics 

An experiment was conducted to establish the best or optimum charac- 
teristics of ERTS black-and-white positives making additive color Images 
for geologic Interpretation, Described herein are both the techniques 
used to establish the optimum Image characteristics and the subsequent • 
geologic Interpretation made therefrom. 

Optimum Photographic Characteristics of ERTS for Geology 

The best additive color presentation of ERTS images was established 
by conducting an experiment In which five different exposure vs. density 
relationships (gamma y) were evaluated. A 5 June 1973 ERTS MSS of the 
WIlIcox Playa area was chosen. Five sets of positives were made from 
this one set of negatives. The photographic parameters of the positives 
were varied to produce a set of transparencies in which the AD (density 
range) of the area of geologic Interest were expanded to as a great a 
degree as possible. Photographica I ly, these were steps 4 through 13 
on the ERTS grey scale at the bottom of the Image. Exposure was care- 
fully controlled in order to keep the geologic detail of the area of 
interest compatible with the requirements of additive color projection. 

Density measurements were made of the negative, release positive 
(y|), and the positive with Increasing gamma. The results are tabu- 
lated In Table 15, in terms of AD of the area of Interest, and the 
expansion factor. The expansion factor Is defined as the degree of 
expansion of the density range of the areas of Interest. 


"r 




86 



Negative 

(Release Positive (y*!) 

Y= I , 3 Pos I t I ve ( nom I na I ) 
y=l.5 '» ” 

y=l.7 " ” 

y=3.0 ” ’• 


D Expansion 
Factor 


Total AD 



Table 15 . Photographic characteristics of ERTS positives reproduced to 
give optimum additive color picture for geologic Interpretation, 


It is Interesting to note that while the total AD follows an orderly 


progression, the D of the areas of Interest, and the expansion factor 

I 

are lower In the y=l.3 positive than In the y=.I release positive. This 
Is due to a slightly different slope of the two curves. The character- 
istic curves that Is negative grey scale density vs. positive grey scale 
density Is shown In Figure 56. 

Screen reproductions of the five additive color Images are shown In 
Figures 57 and 58 . An examination of these images shows an Increase 
In color differentiation with Increasing gamma, with the most dramatic 
Improvement being the Y=3.0, This Is due to a combination of both the 
expansion of the Ad, and the lower Ad min allowed by this particular 
film-developer combination. An examination of the step wedge In the 
grdy scale data block shows the decrease in the brightness Image that 
can be accommodated with Increasing gamma. 

Visual Analysis of ERTS Image for Geologic Information 


A study was performed, employing visual analysis only, to Inves- 
tigate the application of ERTS additive color Images to geology and 
related dJ sel pi Ines, The Images used were made from opt Ima Tl y processed 











Release Positive 
Nom i na I Gamma - 
Nominal Gamma - 


Nominal Gamma 
Nominal Gamma 


Step Wedge Density 

Figure 56. Photographic characteristics of positives made from NASA 
Reiease Negatives to produce optimum additive color image for geologic 
interpretation. 
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Gamma one point seven (y*I.7) 


Gamma one point five (y»I.5) 


Variation of geologic Image quality with Increasing gamma 
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Figure 58 . ERTS image produced using optimum reprocessinq of biack 
and-white positive images. 



gamma 3.0 positives. The two basic types of geologic Information for 
which the Image was examined were; (I) the detection of structural 
llneations, and (2) the dlscriminatlonof lithologic units. 

Geologic, vegetation, and soil maps were visually compared with 
the ERTS image for any possible correlations. Basaltic rock units 

which were recorded on the state geologic map of Arizona showed up 

. ' 

clearly on ERTS-I additive color Images. The Intricate geology In 
tile mountain ranges was not readily detectable. The most recent 
fluvial deposits which were not recorded on the state geologic map 
were imaged and easily seen on ERTS-I. 

A portion of the ERTS-I Image was enlarged and visually compared 
with a large scale USGS geologic map. Contacts between general bed- 
rock and alluvium, and between limestone and slltstone formations were 
obiservable on this enlargement. No faults and other lithologic con-^ 
tacts, present in the study area, were able to be detected from the 
enlargement. 

Objectives of the Geologic Investigation 

The purpose was to determine what useful geological Information 
could be extracted from the additive color rendition of ERTS-I frame 
number E-13 1 7- 1 7222. The two basic types of geological information for 
which the Image was examined were; (I) the detection of structural 
I llneations (faults, fractures, etc.), and (2) the discrimination of 
lithologic units. Due to the scale and resolution Inherent In ERTS-I 

I 

Imagery, lithological discrimination was to be solely based upon color 
differences on the Image and not upon other pa^rameters such as texture, 
topographic expression, and drainage patterns which commonly are 
used in aerial photographic Interpretation. 

^ ^ 


This geologic Inforinatlon can be put Into two separate categories; 

i 

first. Information which has already been recorded on existing sources, 
I.e. geologic maps; and second. Information which has not been recorded 
either because this Information was undetected, or It came Into exist- 
ence subsequent to the published sources. 

Procedure: 

To see how much of the previously recorded Information could bo 
detected on the Image, visual comparison was made between the Image 
and various published maps. If there appeared to be a significant 
amount of Information common to both data, the Information contained 
oJ the published map was drawn on a clear vinyl overlay which was then 
placed over the ERTS Image for a more detailed visual analysis. The 
additional Information which was on the Image but not on the geolo- 
gical maps was drawn o^ separate overlays. The Information on these 
overlays was then compared to non-geo log I cal data, such as, soil and 
vegetation maps. 

Sources Used: 

The main source of geological Information used In this study was 
the 1969 state geological map of Arizona published by the Arizona Bureau 
of I Mines and the U.S, Geological Survey, It would have been most de- 
sirable to have compared the Images with moref^detalled geologic maps 
having a larger scale. Unfortunately, virtually no U.S. Geologic Quad- 
rangle Maps exist for this area. There were a few USGS Miscellaneous 
Geological Investigation Maps In this area but these were poor quality 
and contained only geological ly monotonous areas. The one exception 
was a USGS M,G,|, Map of the southern portion of the Mu|e Mountains, 
Correlations between this map and ERTS Images are discussed In detail. 
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One coincidental advantage if) using the state geologic map was that 
visual comparison with the ERTS image was facilitated because both 
items were at the same scale of 1:500,000, 

Non-geo logical sources included a small scale, very generalized 
vegetation map of the state of Arizona published by the University of 

j 

Arizona, several soil maps which were included In the Soil Survey of 
the Sulphur Spring Valley Area, 1954, and the General Soi i Map of 
Cochise County, Arizona, published by the Soil Conservation Service 
in 1971. 


Discussion of Correlations and other findings in ERTS-I Frame 
Number E-l 317-1 7222 

Since there seemed to be some correlation between the additive 
color image and the state geologic map, all the lithologic and struc- 
tural information contained on the map was transferred to an overlay 
which was then placed over the image; this if? shown in Figure 59, The 
most noticeable of these correlations was the delineation of the two 
major basaltic occurences in the area. These volcanic rock units can 

r. 

be clearly distinguished as having a darker visual signature than the 
surrounding areas. Even the small outlines of these units can be 
easily seen. 

It was found that a large percent of the intricate geology In the 


mountain ranges cannot be detected from the image. The most obvious 


reason for this Is that most areas in the ranges are heavily obscured 


by ve,yet4!tion. However, on close examination, some lithologic contacts 
In the ranges were detectable. For example, the major contact between 
the granite and the schist (Tkg and PCsc, respectively on the geologic 


overlay) in the Little Dragoon Mounta I nsy could be observed on the 
image. No faults could be detected in the ranges. 
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Figure 59. ERTS-I images of south eastern Arizona on which is super 
imposed an overiay of known geology from existing geologic maps. 






Among the additional Information contained on the image but not 
on ths geologic state map of Arizona were the various color subdivisions 
of the basins (Figure 60 ). These different colors could not be checked 
for the possible correlation with different alluvial types because the 
alluvial deposits on the state geologic map were so greatly generalized 
that any such correlation was impossible to make. The vegetation and 
soli maps were also examined for possible correlations. The vegetation 
map was too generalized and the soil maps were too- large-scaled and 
detailed for the resolution of the Image to be of any use. No cor- 
relations were made with these materials. It Is still not known what 
the different colors In the basins, which appear In the ERTS-l Image, 
represent. 

Other features Imaged, but not appearing In the state geologic map 
are the most recent fluvial deposits or alluvium. These recent deposits 
have a white visual signature on the additive color Image. In this 
ERTS freme, they can be seen along the San Simon River, San Pedro River 
and around the WlMcox Playa In Arizona, 

These locations of the most recent fluvial deposits can be easily 
explained. During pluvial periods (times when the area's cl Imate was 
more humid than today) stream water was high and deposits were laid 
down along the major drainage lines. Also, during these wet times, the 
Willcox Playa was a fresh water lake with several streams flowing Into 
It. As these streams entered the standing lake water, they lost some 
velocity and had to deposit their load. Thus, fluvial deposits became . 
concentrated around the perimeter of the Willcox Playa, During Inter-- 
fluvial periods, such as the present, the lake water dries up and evapo- 
rates (gypsum, rock salt, etc.) are precipitated and deposited in the 
center of the playa* 
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Figure 60. Additional geologic information contained in ERTS-I image 
but not appearing on state geoiogic maps. 
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Discussion of Corral at Iqos and Other Findings In the Mule Mountains: 

As previously mentioned. Geologic Quadrangle Maps and MIscella-, 
neous Geologic Investigation Maps were almost totally non-existent for 
this area. The one exception was the Miscellaneous Geologic Invesflgatloh 
Map No. 418, Geologic Map of the Southern Part of the Mule Mountains, 
by Hayes and Landis, published In 1964, 

Due to the availability of this map, a detailed visual comparison 
of this map with the ERTS Imagery was performed. To aid this analysis; 
first, the portion of the ERTS Image which corresponded to the mapped 
area was enlarged to the approximate scale of 1:125,000; then, all the 
lithological Information was drawn onto an overlay, (Faults were not 
Included for the purpose of keeping the overlay uncluttered.) This 
overlay was photographically reduced to the same scale as the enlarge- 
ment and subsequently placed over the enlargement for comparison. 

Again, the two types of geologic Information which were Intended 
to be extracted from the Images’ v/ere the dlscrlrnlnatlon of 1 1 tho logic units 
based on color differences in the Image and the detection of structural 
llneatiphs (faults, fractures, etc,). 

Lithologic discrimination of this area was greatly limited to the 
observable contacts between the alluvium and the general bedrocks of 
the Mule Mountains, This all uv I um/bedrock contact Is clearly visible 
along the total eastern side of the Mule Mountains (see ( I ) on the over- 
lay)/ and also In the area SE of the major ta 1 1 Ings pond, (2) on the 
overlay* Only In one place Is the contact between the valley fill and 
the pediment gravels (Qts and Qg on the geologic map, respectively) 
observable. This area Is east of the San Pedro River where the Green- 
bush Draw Joins It (3) on the overlay). However, going away from this 
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LITHOLOGIC OUTLIN'E 
OF Tf' i>Ol. irir,«L' P/- 


THE G.,:-Lv)r,ic p .p . 

THE MULE MTL. . ^ 

Figure 61, Lithologic outlines for the Geologic map of the Southern Part of the Mule Mountains 
Arizona, superimposed on an enlargement of ERTS additive color rendition of that area. 




area, the contact becomes completely unobservable. 

’ t 

The Juniper Granite (J| on the geologic map) can be seen on the 
enlargement (4) but this Is because of the growth of the vegetation 
around, but not on, this lithologic unit. The outline of the Juniper 
Granite also can be clearly traced on the 1958 BIsbee 15* topographic 
quadrangle due to this same preferred growth of the vegetation. At 
two separate locations In the southeastern portions of the Mule Moun- 
tains, contacts between I lmestor>9 (5) and slltstone formations are 
visible^ as color differences on the ERTS Imagery, the limestone being 

i 

whiter than the slltstone. 

No other bedrock contact or faults and fractures could be detected 

1“ 

i ■ " ' 

on. the enlargement. One possible reason for this might be that vegeta- 
tion covers the areas where most of these contacts and faults occur, 
l,e, the center of the Mule Mpuntalns, 

The most prominent feature seen on the enlargement is the tailing 
pond (6) and the smaller ore dumps and other mining waste sights, located 
just southwest of BIsbee, In every case, these waste sights appear to 
be larger on the Image then on the geological map. This discrepancy 

coiild have different explanations. These waste sights may have been 

I 

expanded since the geologic map was published in 1964, This would demon- 
strate ERTS Imagery as an effective tool for updating geologic maps. 

The only other explanation Is that the Image Is not just showing the 
area which Is physically occupied by the waste sights, but Is showing 
ths total area which Is affected by these features: (Areas where vegetar 
t Ion may have been killed by adverse fumes or by acid seepage Into the 
groundwater) ^ case, ERTS Imagery would be an effective source 
of data for environmental Impact studies of mining dumps, 

' 
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The enlarging of the ERTS Image aided to some degree the Interpre- 

' i ’ ■ 

tablllty of such,! The alluvium/bedrock contact SE of the tailing pond, 
(2) on the overlay, was not detected on the original Image, ' However, 
the enlargement did not Improve the detection of many faults and other 
lithologic contacts present In the Mule Mountains, 

Conclusions! and Recommendations: 

Except for the basaltic lava^flows, no one lithologic unit could 
be completely mapped from the ERTS frame E-I3I7-I7222, This lack of 
ability to discriminate between different rock types may be attributed 
to the facts that: (I) the mountain ranges, which contain practically 
all of the outcrops of the area, were heavily covered by vegetation, 
and (2) the bedrock in the basins was covered by thick accumu lations 
of alluvium. Furthermore, the possible correlation between different 
types of a I luv I urns and different colors on the image was made impossible 

X 

by the lack of good supporting materials. It Is therefore suggested 
that, if future studies of this nature are to be conducted, an area 
should be chosen In which there is a minimum of vegetation and super- 
ficial deposits covering the actual rock outcrops, and where there has 

been a considerable amount of detailed geologic mapping. The latter^ 

would assure the availability of ample amounts of good supporting 
materials which, in turn, would reduce the need for field checks. Also 
aerial photographic Interpretation of the area may be used In place 
of field checks. 
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Section 6 ' 

CONCLUSIONS 

Image Characteristics Related to Reflectance Spectra 

The image density in the ERTS-I release positives was found to 
be indicative of the percent reflectance of the Wiilcox Ptaya, Arizona, 
Percent reflectance data of the twelve sampling stations within the 
Playa were averaged and compared to the mean value of 250,000 digital V 
brightness measurements of the Playa obtained from ERTS-i release 
positives. These data are shown in Table 16 below. 


Ba nd 

In Si+u Spectra 

ERTS Image 
Brightness 

~ s 

Incident 

yw/cmVnm 

Percent 

Reflectance 

Reflected 

Uw/cm*/nm/sr 

500-600 

79 

32 

25 

124.9 41.2 

600-700 

72 

38 

27 

155.4 45.0 

700-800 

60 

40 

24 

156.1 44.4 

800-1 100 

1 

35 

42 

15 

152.8 34.0 

Table 16. 

In sltii spectr 

•al ref 1 ecta nee measurements and 

correspond i ng 


ERTS-I image density measurements. 

These data were normalized with respect to average values at 950nm, 
The average digital brightness measurements were norma I i zed to Ba nd 7 
and the average percent In situ spectral reflectance data normal I zed to 
950nm, The results are plotted In Figure 62. 







normalized to 950nm and ERTS-I image brightness data normalized to 
Band 7. 

While normalized ERTS-I image brightness is related to spectral ref- 
leLtance measurements normalized with respect to 950nm, calculated radiance 
values (viz, incident irradiance times percent reflectance) showed little 
correlation as can be seen from the above table which- probably due to 
.the effects of the intervening atmosphere, 

A comparison of the spectral reflectance data obtained at the Prescott 
National Forest and the In situ percent reflectance of Pinion Pine, Juniper 
and Manzan I ta showed no Relationship whatever even though these species com- 
pletely cover the gr^^ of the Copper ’Basin Area of the National Forest, 
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Analysis of the ERTS-I Image Brightness 

Statistical analysis of 294,912 digitized pixels of ERTS-I positive 
images showed that: 

- Ihe WII Icox Playa Is darker In Band 4 than any of the other 
ERTS-I bands, 

- the brightness of the Willcox Playa Is approximately the same 
In Bands 5, 6 and 7, 

- the variability of the brightness data is less in Band 7 than 
in any of the other ERTS-I bands. 

- the frequency distribution of brightness has tv/o modes (peaks) 

In Bands 4 and 5 which are absent from Bands 6 and 7. 

I , _ 

Statistical analysis of multidata ERTS-I negatives and positives 

showed variations in the spectral reflectance of the imagery of 
the Willcox Playa between spectral bands and In the same spectral band 

. V 

on different dates. As is discussed below, this difference In Image 
brightness between dates was due to standing water. 

A statistical analysis of the difference In mean digital brightness 
measurements of cloud tops and shadows showed that the image contrast 
in each band increased as a percentage of Band 7 as foUows: 

X 100 = 74^ - - 

Band 7 

X 100 = 19 % 

Band 7 

X 100 = 88^ 

Band 7 

The digital brightness of ERTS-I Images of the Prescott National 
Forest Is not In any way related to the spectra of tree species growing 
there. It was not possible to differentiate Manzanita, Juniper and 

■ ■ ' t . ' 

Pinion Pine. ■ 
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Statlsticah analysis of digitized ERTS-I negatives showed that the 

. * » 

mean brightness difference between Images of terrain In sunlight andjn 


shadow expressed as a 


percent of that difference In Band 7 Is: 


X 1 00 = 34JS 

Band 7 

i' 

X 100 = 60JS 

Band 7 

Q^nd 6 X 1 00 = 16 % 

Band 7 

i ■ 

' I 

It was also shown that these differences Increased with decreasing relative bright 

ness. A nomograph to correct for atmospheric conditions was derived using cloud 

I. ,, ' '• . 

tops, shadows and adjacent illuminated terrain areas. (Figure 54 on Page 85), 

Additive Color Presentation and Analysis of ERTS Multidate Images 

The dynamics of spectral reflectance changes over time were shown 
toj be capable of presentation on an additive color viewer. The purpose 
was to establish if differences in standing water in the Willcox Playa 
could be detected in ERTS-I imagery taken on different dates. 

ERTS-I images of the Willcox Playa acquired oi, September 27, 

October 14 and November 2, 1972 were used. The image of the Playa 
was located in different positives of the frame for each date so that 

I 

only the Intersection of the three Images could be combined in color. 

Reprocessed positives were made from the NASA release negatives. The 
September 27 Band 5 was projected as blue, the October 17 Band 5 as 
green and the November 2 Band 5 as red. A copy of the recombined 
Is shown In Figure 63. 

Ordinarily, the colors in a mul t ispectral display are indicative 
of reflectance differences which exist as a function of wavelength. 

However, the vivid color differences which are seen In the multidate 
ERTS Image bn the next page, represent; reflectance changes In the 
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Figure 63. Multidate additive color image of the Willcox Playa 
showing changes in standing water, 

600-700nm spectral band as a function of time. The colors within 
and surrounding the Playa are Interpreted as follows: 

- A colorless shade of white or grey indicates absence of any 
standing water between 27 September and 2 November 1972, 

- Yellow indicates areas where a decrease in standing water 
occurred between 27 September and 14 October and little change 
from 14 October to 2 November. 

- Orange also indicates a decrease in the amount of standing 
water between 27 September and 14 October, coupled with a re- 
duced rate of decrease to 2 November. 
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- Cyan are areas where standing watpr changed little between 

27 September and 14 October, then Increased greatly by \ 

2 November. 

-Red Indicates areas where surface water did not change between 
27 September and 14 October, then decreased between 14 October 
and 2 November, 

- Green Indicates a reduction In surface water from 27 September 
to 14 October and an Increase between 14 October and 2 November. 

- Blue Indicates the opposite effect from those areas that appear 
yellow; I.e,, the amount of surface water Increases between 

27 September and 14 October with no change between 14 October 


and 2 November, 

These procedures for using additive color analysis of multidate 
ERTS Images developed under this study was the Interpretation technique 
used to map stand i ng. water during the 1973 Mississippi Flood, A 
summary report of this flood mapping project using mult I date ERTS 
Imagery appears In Appendix G , 
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Section 7 


New Technology 


The research reported herein Involves the mu I tl spectra I analysis of 
ERTS scanner data. 

The significant theories and the description of the Interpretation 
techniques used have been given In Sections 2, 4, and 5 of this report. 


Appendix a 


"SYSTEM 800" 
DESCRIPTIVE LITERATURE 


A typical System 800 interactive display system using the Model 805 Color Display combined with the- 
Model 810 Digital Processor is shown on the opposite-page. 

t 

COLOR DISPLAY , 

The large screen color picture display is constantly refreshed at 60 fields per second from the refresh store 
and is constructed from a rectangular array of 512 by 512 picture elements. Each element, under individual 
program control, may be displayed in 32 colors, in a 32 level grey scale, in black, or in white. ... 


COLOR KEYBOARD 

The push-button keyboard in front of the display provides operator control of the picture elements accord- 
ing to their value. Thirty-two columns of keys, one column for each color or grey scale value, change the’ 
disf)lay of all picture elements of a particular value from color to grey scale or to black. 

In addition, under keyboard control, one or more equal-value sets of elements may be combined to reduce 
the number of colors or grey shades in the display. 

Thus, the Color Keyboard speeds the analysis of displayed data by providing immediate operator cOntrbV 
over the number of colors and grey-scale contrast. 


DIGITAL PROCESSOR 


The Processor and Magnetic Tape Units are supplied to meet the needs of specific applications or 
805 may be supplied alone with the proper interface for the customer's processor. 


the 


Model 


Operator control is provided through the Teletypewriter keyboard and results and programs are listed on 
the teleprinter. Program entry is rapid and convenient using the high-speed Photoelectric Tape Reader. 


REFRESH STORE 

The refresh store contains the 262,144 picture elements and continuously transmits them to the display at 
a peak rate of 10 million per second. At the request of the Processor program, the picture elements are- 
transferred directly to or from core storage in blocks of 512 picture elements (one horizontal line). The 
transfers take place over the Processor input-output channel in a similar manner as other input-output 
devices, the magnetic and paper tape readers. ; - ' ’ , ’ 

OPERATING SPEED „ . , - . 

When operated with the System 800 Digital Processor, the entire 805 Display may be changed- in less than 
eight. seconds, new picture data may be added at the top or bottom and the picture moved up or down 60 
horizontal lines per second, or the Display may be changed from a positive to a negative, from color to 
grey scale, or from picture to graphics only, in 1/30 second. The average access time for 512 elements is 
1/60 second. 

















SYSTEM OPERATION 


Tlie block diagram on the opposite page shows the major: components of the complete Color Display 
System. The Model 805 Display is shown to the right of the input-output (I/O) bus of the Model 810 
Processor. 

Tl|ie high resolution color picture refreshed at standard teievision rates is constructed from a rectanguiar 
array of 512x512 picture elements. Each picture element is specified by a 5-bit picture code and a 2*bit 
graphics code. The 32 picture values and the four graphics values are individually controlled from the 
processor to provide a wide variety of displays. 

Two data paths are shown to the 805: one through the upper interface electronics to the Refresh Store for 
transfer of the picture and graphics codes, and one through.the lower interface electronics to the Display 
Control for transfer of control words, status words, and cursor address. 


Tfie digital-to-analog converters in the Color Generator convert, under control of the Color Keyboard, the 
picture and graphic codes, which are continuously being circulated in the Refresh Store, to red, green, and 
blue video signals for the Color Picture Tube. According to the four possible values of the graphics code, 
the 32 possible values produce, for each picture element, either one of 32 colors, or one of 32 grey shades, 
or black, or white. The picture and graphics codes are separately accessed through the Buffer Store so that 
they may be modified independently. 


piGITAL PROCESSOR > ' 

The Model 810 Digital Processor and Model 810-1 Magnetic Tape Unit are designed to operate with the 
Vjodel 805 Display as a stand-alone system for the analysis of picture information. The system can perform 
data editing and formatting, data analysis, and the production of graphics. Data is normally first entered by 
the 810-1 Magnetic Tape Unit. This data may have originally been digitized from a sensor (such as a ther- 
mal, radar, or sonar scanner) or may have been generated by a large-scale digital computer or the Model 810 
Itself. 

j 

The data on the magnetic tape is read, a block at a time, into the processor core storage, and, if not already 
formatted in a 512x512x5-bit matrix, is converted to 512x5-bit blocks of horizontal rows of picture- 
elements for output to the Display. The data is transferred to the Buffer Store and then to the picture code 
section of the Refresh Store. The latter transfer takes place by a software-controlled output to the 805 
control of a control word which contains the address of the line in the Refresh Store, 

7Jhe transfer from the Buffer Store is made to one of the 512 horizontal rows of the Refresh Store specified 
bjy the control word placed in the Display Control by a programmed transfer. 

The graphics codes, 2 bits for each picture element, are transferred to the Display in the same way as the 
picture data. Areas may thus be shown In grey-scale or color depending on parameters calculated by the . 
Processor. v 
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REFRESH STORE AND BUFFER 

The picture and graphics are stored digitally in the Refresh Store. The horizontal lines, each containing 512 
picture elements, are circulated, even lines followed by odd lines, so that they may be continuously read in 
synchronism with the standard interlaced scanning pattern of the color picture display. Successive lines (Y 
addresses) are read every 1/60 second. * . 

The Refresh Store is addressed by line |Y) address from the control word and, at the proper time, the com- 
plete line of 512 picture elements is transferred to/from the Buffer Store. The Buffer Store then transfers 
to/from the Processor core using a high speed block transfer. 

i ■ ; , 

The Buffer Store may be used for transfer of picture or graphics data as specified by the control word sent 
to the Display Control. 


COI.OR GENERATION 

The 5-bit picture codes and 2-bit graphic codes are read continuously from the Refresh Store. The Color 
Generator using special digital-to-analog matrix conversion circuits produces red, green, and blue color 
video signals to drive the three electron guns of the Color Picture Tube. Colors are produced by generating 
proper proportions of the red, green, and blue primaries on the screen. The color produced for each picture 
element is controlled by the 2-bit graphic codes as foliovys: 

Binary 00: One of the 32 standard colors is assigned according to the 5-bit picture code (color 

picture). 

Binary 01: White is generated regardlass of picture code values (v>»hite graphics). 

Binary 10: ’ Black is generated regardlt?ss of picture code values (black graphics). 

Binary 11 : One of 32 shades of grey ranging from black to white are assigned according to the 5- 

bit picture code (B/W picture). . . 

Thus, the graphic code determines the type of display for each element. 


COLOR KEYBOARD 

In addition to the software control of the display the colors are also under control of the Color Keyboard. 
Four ipterlocking push-keys are provided for each of the 32 colors and control the display of each color as 
follows; 

■ i . 

1) Normal: Normal Operation, the picture-elements of this value are displayed on the screen as 

! determined by the picture and graphics codes, ' . 

2) Black: Every picture-element of this value is now displayed in black. 

3) B/W: Every picture elementof this value is displayed in grey scale. 

4) OUT: Every picture element of this value is displayed in the same mode as next higher order 
value. For example, the OUT key of value 10 displays value 11 in place of value 10 so that 10 

I and 1 1 are now the same color or grey level. 

Th^ black and B/W keys perform the sarhe functions as the graphic codes on ^ elements of the same value 
or color and are a great help to the operator for recognizing particular colors and features. ; 

The OUT keys eliminate colors so that the complexity of the color display may be reduced quickly without 
need for the processor to rewrite the picture data. 
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SYSTEM CONTROL 


The system control 'accepts a control word from the Processor and supplies a status word to the Processor, 
When the cursor and joystick option is provided, a two-word X-Y address is also supplied to the Processor 
indicating the position of the cursor. ' 

Program control of all picture elements simultapeously by a single 16-bit control word produces an all color 
display, an all grey scale display, graphics only, a complemented (negative) grey scale display, and roll up 
or. roll down of entire picture. The latter command provides a rolling raster for display of continuous strip 
pictures. 

The control word consists of an operation code and a line (Y) address. Only the Buffer to Refresh Store 
transfers require a line address. Specified bits in the operation code determine whether: 

1) elements are displayed as specified by the graphics codes as described above. 

2) all elements are displayed as specified by the graphics codes except color elements are displayed 
in grey scale. 

3) all elements are displayed in color (without graphics control). 

4) ail elements are displayed in grey scale (without graphics control). 

5) the graphics black and white elements only are displayed (without picture data), 

I 6) any of the above are displayed with the picture data complemented. 

i 

These control functions require only one I/O instruction for operation and immediately switch the display 
from color, graphics, to grey scale under program control. Thus, the program! may present ten types of 
information display to the op.erator without reloading the picture store. 

In addition to the above operations, the operation code includes roll-up and roll-down. These operations 
immediately move the displayed picture up or down two horizontal lines. Continuous movement of the 
picture may be performed by writing new data at the top or bottom of the picture before each roll com- 
mand. In this way new picture information can continuously replace the old, producing the effect of a 
continuous roll. 

t 

The status word provides the Processor program with updated information on the display status. An inter- 
rupt is given to the Processor at the completion of each operation that requires more than one Processor 
instruction time to complete, such as a buffer-disc transfer. 

■\ 

i • ' ' 
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SOFTWARE 

The Mpdel 805 Display is provided with software on punched paper tape for operation with the Model 810 
or customer's Processor. Included with the software package are: 

1) Diagnostics for the Processor, core memory, and display system 

2) Machine language assembler 

3) Basic compiler 

4) Load and Utility routines 

5) Debug Program 

6) Math Functions • • 

The Processor is supplied with a software library including'display diagnostic routines for checking the 
operation of the processor and display unit. User programs may be entered originally, edited, and assembled 
using the teletype keyboard. Programs are stored on paper tape and quickly loaded using the high-speed 
photoelectric tape reader. Operational activities are controlled through the teletypewriter keyboard and the 
Processor front-panel switch register. 


SPECIFICATIONS 




DISPLAY 

TYPE: 19" Shadow Mask Picture Tube, 670x500 color-dot trios. 


\ 


COLOR DISPLAY: Number of Colors - Selectable 1 thru 32. ^ 

Color Assignment - Any one of 32 colors assigned to any picture value. The colors are four shades of 
Yellow, Cyan, Green, Orange, Magenta, Violet, Red, and Blue. Colors are controlled by graphics data 
and keyboard to produce black, white, or monochrome. • 


MONOCHROME DISPLAY: Number of Grey Levels • 32. Monochrome superimposed on color under key* 
board and program control. 

graphics DISPLA Y: Graphics displayed in black or white, Graphics display Cinder program control. 
/?£SOZ.L/77G/V; 512x512 raster points stored, 512x489 displayed. ’ 

6E0/WfTR/CG/ST0/7r/0/V; Less than 2% overall. 

S/Zf; 16x12 (4/3 aspect ratio), 12x12 (1/1 aspect ratio). 

REFRESH RATE: Standard 525-line television interlaced 2:1 , 60 fields/sec,, 30 frames/second. 

DISPLAY KEYBOARD: 4x32 pushbuttons for selection of color, black, B/W, and Out. 

COLOR PATCH PLUG: 64 pin plug for change in color order. Prewired plugs for normal and reverse 
order supplied. „ 

CURSOR: Joystick control of cursor position. Activate switch. Maximum interrogation rate 30 times/ 
. second. 

MODEL B-IO DIGITAL PROCESSOR 

CORE STORAGE: 8192 16-bit words. _ 

NUMBER OF REGISTERS: 8 

INSTRUCTION SET: Over 400 hardware instructions. . 

SPEED: Memory Cycle - 1.2 microseconds, Add/Subtract - 5.2 microseconds. 

PERIPHERAL EQUIPMENT: Teletypewriter ASR 33, Photoelectric Paper Tape Reader (300 chars/ 
second). Model 810-1 Magnetic Tape Unit - 9 track’, IBM compatible, 800 bpi, 37% ips, lOVi" reels. 

SOF7Wi4/?£’:Diagnostics,MathFunctions,Assemblers,BASIC,Debug.- 

■ ■ E ' , - V ■ " 

POWER AND PHYSICAL 

SYSTEM POWER REQUIREMENTS: 105/125 volts, 60 Hz, 20 amp circuit. ’ 

£/VWflOA/Aff/Vr; -H5° to -i-30‘’C, 10% to 90% Humidity* 
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• RBN 


i'^RUN C-rSPLV 
li i #P 2 


$RUN MARKER " 
i $RUN HISTOR 
[SCAN : 

i Z = 202 AT <283, 438) ' 

I Z=190 AT <283. 448) 
iZ=t99 AT <284,428) 

I Z = 204 AT <276, 428) 
IZ-204 AT <288, 440) - 


ZOrjA 


ioori;-! 

OP Pt^AVA , 
Bi^rvyO ^ ^ 

18 rim 










I >RU« flREfl 

1 coroi^ 0 
I GRftV; ' 032 


yS PEk. )Vl. peU£>^rm 


2 3 

096 ‘ 128 


PERCEMT: 28. 92 24. 40 24. 72 15. 26 10. 72 03. 60 00. 28 00. 00 


f RDW DISPLV 


•RON DISPLV 


■■--r&n c § 







APPENDIX D 

"SYSTEM 800" 
COMPUTER PRINTOUTS 


PHOENIX, ARIZONA 
SPECTRAL ANAlYS IS 





$RUH DI5PLV 


ftlPO - 


tRou "Drcrir 


S 00 vs»rB/^ -.comPu.r.C'/z., 

Pfeiivroo r ^ V I 


$RUN RDDER 
HP1VP0 ' ' 


•R0N DISPLV 


JP. . _A 0 ni a_„ 

A iJ At-v| ^ 1 6 - 


^RUN 

itPO 


!|R0K negrte 


!$r|u» 

]»P0 


DI5PLV 


• RUM D'l SPUV' 


.E£3^.-i Po^iTj/i O P tifiobi 




•R0N DI5PLV 


I RUM 


IR0M 


MARKER' 

OIQIT 


$RUN HISTOR 


SCAN 


$RUH AREA 


0-„ TB^^ ij. A ^ kjr\)k^to OiJjjTi 2. 

„ „ „ , J ^ L4 i hJ T l)% iArw E. 



frlk}^. /V ^ . -- 


COLOR: 

GRflVl^ 

EflfCEN 


14. 7 


0^ ■V<«=^1'28 .^1^0 

2. 24^7. 7^8^" 29. 05. 42^00..‘3l 00.*b0 ... 00. 00 4-.-— 







. KI 

fRUir"MFIF:KER- 
fRB^^ [>IGIT~ 

I 

♦ R'tJff HISTOR 


r 

Z*222 

Zs202 

z=2ie 

Z=204 
Z = 265 
Z=20e 


fir <168, 

AT (169, 
AT (IZl, 
AT <1Z5, 
AT <162, 
AT"<162, 


292> 

300> 


250> 

244> 


Z = 0Z? 

AT 

<152, 

24 6 > 

Z = 069 

AT 

<153, 

238) 

Z=,055 

AT 

<151, 

248) 

Z='663 

AT 

<151, 

24 2) 

Z = 059 

'AT 

<152, 

238) 

Z=0?9 

AT 

<155, 

232) 

Zs0G8 

AT 

<155, 

282) 

Z = 084 

AT 

<155, 

290) 

Z = 0Z1 

AT 

<154, 

284) 


J?/i 5 /<vC> b. 

im. 


CLoV'(^ _SAa;^ 




Pe^S'jrMeA 


COLOR: _ 0 1 ^ 2 3 4 5 6 ? 

GRftV: 032 ■ G64 096 128 160 192 224 25 

PERCENT: 01.98 04. 54 12. 62 24.80 31.82 20.60 03.50 O0. 08 


• RBN C-rSPLV 


^ “ t fl 

Prf^/liL 1^1 Z 

Phhitj I X 2-0 A 


Olo\)^ pA-rvf^ 7 

























APPENDIX E 


"SYSTEM 800" 
COMPUTER PRINTOUTS 


PRESCOTT, ARIZONA 
SPECTRAL ANALYSIS 


INV CUD! 

$ 


P/h ^rooT 


RUN t>ISPLV“ 
#P2/S 


<fpUN DISPLV 
»P3 


_S|P S^’Vyt*e^ A-N i^u. ..OP 


.HoPPER. A Ai.i/0j P^BLt^TT Ai^iBo uA 


$RUN HISrOR'- 
SCRN 


1 leuNcnsER 


1 

PliBi.^TT ] 

1 IMf CMD! 

^ ■ fl 

2 9 {^UL 1^1 Z' iv^a/L Ti o/J 

1 * 

■ i 

1 RUN RRER 

’ ' ' 

I ' ' ' ■ ' " ‘ 



COLOR 

or; 

Trceh 

!, 

$R!UH 
5 CAM 

z-oee fl 
couHr= 

COUNT= 
CjOUNT = 

' ’couNr= 
couNr= 

0 > COUN 
CpUNf= 
■ COUNT = 
J CpUNT= 
^~COUNT = 
1 > COUN 
‘CpUNT = 
COUNr= 
' cbUNT= 
2 >' COUN 

■""couMr= 

COUNT= 

~couNr=: 

COUNTS 
3 > COUN 
COUNTS 
“COUNTS 
COUNT= 

COUNTS 

COUNTS 

COUNTS 


121 

RT 

RT 

RT 

RT 

RT 


:T0R 


T <122 
002168 
002621 
002282 
001422 
001166 
T= 001166 
002422 RT 
002821 RT 
002280 RT 
002161 RT 
T= 002161 
002469 RT 
002861 RT 
0021Z6 RT 
T= 002176 
“0015 78’ RT 
001162 RT 
001229 RT 
001526 RT 
T= 001526 
000820 RT 
000650 
000449 
000270 
000264 
000270 


049 

058 

047 

027 

022 


RT 

RT 

RT 

RT 

RT 


GRRV. 
GRRV. ■" 
GRRV. 
GRRV:*~ 
GRRV. 

RT 022 GRRV. 
054 GRRV. 
060 GRRV; ~ 
062 GRRV. 
064 GRRV: 

RT 064 GRRV. 
082 GRRV. ■ 
091 GRRV. 
096 GRRV: — 
RT 096 GRRV. 
118 GRRV7~ 

127 GRRV. 
120 GRAV.*'" 

128 GRRV. 

RT 128 GRRV. 

142 GRRV. 


14 A_S. tr.. w£._.OE_llinL 
J rJm. ^ TO 

PlUUliAf A. 




152 

158 
161 

159 
161 


GRRV7 

GRRV. 

GRRV.' 

GRRV. 

GRRV. 






^"cWnT="" 

iiW370" 

AT’ 

Ibl 

"GRiRVT 

COUNTS 

000367 

RT 

162 

GRRV. 

" *COUNT = 

000364" 

RT 

159 

GRRV: 

COUNTS 

000370 

AT 

161 

GRRV. 

" COUNTS 

000576 

AT 

160 

GRRV." 






4> COUNT= 000576 
' C0UMT= ' 000183 AT 
C0UNT= 000059 AT 
■' C0LIMT= 000182 AT 
5> C0LINT= 000182 
C0UNT= 000082 AT 
000087 
000056 
000034 
000056 
000056 
000088 
6> CCiLINT= 000088 
' C0UNT= 000024 AT 
C0UNT= 003002 AT 
?> COUNT = 003002 


AT 160 GRRV. 
184 GRAV: 

191 GRRV. 

192 GRRV: 

RT 192 GRRV. 

209 GRRV. 


COUHT= 
COUNT= 
C0LINT = 
COUNT= 
COUNT- 
COUNT= 


AT 

RT 

RT 

AT 

RT 

AT 


216 

221 

223 


225 


GRRV. 
GRRV.""'" 
GRRV. 
GRRv: "■ 
GRRV. 
224 GRRV, 

RT 224 GRRV. 

254 GRRV: 

255 QRflV. _ 
RT 


255 GRRV. 





! IfRUN 

RRER 

■ 


7 



i V :? 

1| COLOR 
1 iGRflV: 

0 

032 

1 

0 64 

2 3 

096 128 

4 5 6 

160 192 224 

7 

255 


r iPERCENT; 07.00 

1 ! 

‘1 i 

1 ■ i 

23. 38 

36. 10 20. 06 

08. 56 ’ 02. 5'2 00. 84 

01, 44 

* 


B ^^RUN 

1 ;■ 1 

MARKER 






f "!tRUN 
i #P0 

OIGir" 


• . . 




ij ;rRUN 

HISTOR 


• 




:r 'SCAN' 




. 


'f 

1 "IRUN 

i ;C: 

RRER 



• 



i i -.1 




P&UBfjr MeA 

^ ' • ■■■ 


1 ^ OLOR 

:: 0 

1 

2 3 

4 5 6 

7 

’■ ■■■ '4 '■ 

r GRRV; 032 

1 PERCENT: 04.34 

064 

11. 96 

096 128 " 

27. 54 29. 34 

160 192 "224 

15. 50 07. 18 03, 16 

255 
00. 90 

'a 

■■ . 

|l $RUN 
ft-ttpe— 

y 1 ■ 1 

li ^RUN 

OrSPLV 





' j 






,;.V' 

Si 

li 1 

MARKER" 





■■■ 'i'i' 

r iRiiiN 

r j 

DIGIT' 






I .rteLiN 

HISTOR 




t ■ ■■ ■ . .. 

■ ■ ■ . -p' 

f sfcflH" 

iPRUN 

RREfr — — " 

~:;r - [ 




■’.'i 




-COtORT 

GRRV: 


0 - 


2 — 2— "6 7 

064 096 128 160 192 224 255 

PERCENT; 09. 64 16. 80 25. 84 19. 28 11. 82 06‘ 92 ' 04. 00 " 05. 60" 


032 


•RUN OISPLV 








sar.:-: 






'll 


IRON NEGATE 




=«»5«se^'aKie^. 




.« 


• 


$RUN AODER 
t » P 2. P 3 


IRON DISPLV 


ipN NEGATE 

I 

IRONPfiDDER 


8 qO ( IcJnnfyM .^ 

^JTbg C-f 


IRON DISPLV ’ 


$ 


RUN HI STOP 
SCAN 


S > Pi 2 1 TitAu, Av. A<^u s A„_a/=l 


A TmK>LA:l±t^ ^ E fi/r£(J%. 


rRUN^MARKER 
i;RUN HISTOR 

Scan 


t , 


Pr 


fRUN MARKER 
$RUN HISTOR ” 


SCAN 


- I ,pO£.m\>g. > rnvVO C.C ovjD ^ ^ 

Ai^o .. „„ U ;> LUC / Q ^ A - y ^ r 2 -DyA 


!> V . 


•iRUN MARKER 
^RLIN HISTOR 
SCAN 


4 / i£Pr. 1971 ; /Co(a C)H>ri2/5TLo^ 


$RUN HISTOR 
SCAN 


0 TlWv li A S.ET uP TP bt (a»r/Z t^ /^ML 



2=220 
2 = 2 26 ^AT 
"if = 22 /' AT 


C o\jr 0-^ ^ or>f5^ 

/WCi TO 


2=^21 AT 


^"229 AT 


Ti 4 ^r Al>l 6/ v ^' D ^ FAi-u wi 57 H*wL 

■ • 4 .'. 

O - ZS^ O' t> I TlTATV)\i i.’cj^^ X. 


La : 




5^ 


■ ■•'t'" 

‘ i» 



jj rj r-j rj rj r-j w ro r-j 




f 



t/i6 Ui .JL^W^Ui i,o /Vt^u /^/WOl "• 

rVvu. Orzsr^ :bri^iiAnt//U. 

, uiVu^ ... ■ ■ ..- 


SRUM MRRKER 
rRWN*HISTOR 
5CRN 


•i ~ 







b f b )P ^4} \JM^ , fiO^ — - 

^ 'Fbvj^ /E^LTI-I iS/WbSk 

LP’^kCEMT 0 0 01. 38. 5 jt.^-^-c--0'0^"9 S'^^.».«..-.^^'f 


^RLIH MRRKER 
iJR^UH-HISrOR 


. £^Pi. "•! 

yiu>r;/975^ 

tAJ^LLrl/)iL 


Z = 069 
I Zf oee 
i'z=05e 
I z«05r 
rz^040 
[_Zs056 
I rs056“ 
! Z«D44 


(117, 

(117, 

(119, 

(119, 

(119, 

(119, 

(114, 

<114, 

(1107 

(108, 


260> 
254 > 
24 8 
240) 
234) 


214) 

208) 

222 ) 


/rJ 

|i/^nrs/Jb^ S’ 


\J\^ t. /^ L 

*•■■ 

IO\3 TTi/, rAOWiTofl/ 




,-^2K 


44 





C 0. 


V2 


f" / 


yi 


=.:.0.-r--. vLQ.R-'3.. .?5.R). 


(VUT* / 

f>f>,iK,fc..r'j; 





2^00 
2 = 144 
2 = 145 
2 = 135 


2 = 146 
2 = 143 


- «■ MI$X>V»3ISSl|««Du 


igar.w^ 


<065, 238> 
<085, 230) 
<086, 220) 
<090, 216)‘ 
<090, 212) 
<085,"210>' 
<116, 220) 




.. BAT4b S’... pEfc-tCrOr 


$RUH RREfl. _ 

COLOR; "1 ■ 1 2 

GRflV: 032 064 096 

PERCENt : 01. 56 01. 80 04. 80 


3 4 “s"" -"~~6 7 

128 160 192 224 255 

09. 42 49. 28 *30. 34 02. IT'™^ 075 8" 


RUN' MARKER 
♦RUN HISrOR 
SCAN ” * 


JaJ J ^ ^ Co >t „ 


2 = 229 

AT 

<154, 

228) 

2=233 

AT 

<154, 

220) 

2=231 

AT 

<152, 

260) 

2=231 

AT 

<158^ 

252) 

■2=230" 

“At 

< 156, 

24 0> 

2 = 226 

AT 

<153, 

234) 

2=233 

AT 

<151, 

250) 

2=230 

AT 

<144, 

252) 

2=234 

AT 

<142, 

266) 

2 = 232 

AT 

<141, 

228) 





C. ’ 
•y 


C.- 


ISS'flT” 
143 AT 
145 AT 
147 AT 


Z = 


143 AT 

144 AT 


c. ■ 

Z=143 AT 


f 

f Z- 

z- 


145 AT 
145 AT 
145 AT 


'C 0 9 0, 
<.‘ 088 , 
C085, 
C085, 
C088, 
C087, 
(0 91, 
C091, 
<089. 
<086, 


274) 

264) 

254) 

242) 

232) 

224) 

216) 

208) 

226) 

240) 


f\/o IT 


•W . 


$ RUN AREA 

COLOR: 

GRRVV 


PBMk}T JidEA 


PERCENT: 


0 

-032- 

01. 30 


1 

06C 


2 

"096' 


3 

128 


4 

160 


5 

192 


6 

224“ 


00 46 01. 62 09. 64 62.80 20. 92 02. 34 00. 82 






JRUM HISTOR 


SCAN 

Z = 2.U AT 
2=2^0 AT 
Z-2'iD AT 
Z=2|2 AT 
Z=229_AT 
Z^228 AT' 
Z=22? AT 
Z=226 AT 
2=228 AT 


Tl52‘. 

d52, 

C152. 

<152, 

<152, 

<151, 

"<1447 

<143, 

<138, 

<136, 


2fsy 

268) 

264) 

254) 

248)~ 

240) 

240)" 

254) 

262) 

272) 


* 7 sgpr. 


_0dU:).U.5- 


Z=067 AT 
2=051 AT 
“2 = 052 AT 
2=058 AT 
2=058 AT 
2=042 AT 
2=035 AT 
2=050 AT 
2 = 045 "AT 
21=054 AT 


<113, 

< 112 , 

<117, 

<117, 

<114, 

<114, 

< 112 , 

< 111 , 

<109, 

<107, 


264) 

256) 

250)" 

244) 

240) 

234) 

226)' 

224) 

220 )' 

216) 


ftj 

SI4A^0W___ 

IE ^ ( j 


2=t:45 AT <093,206) 


2=153 

“2 SI 52" 

2 = 156 
2=155 
2 = 156 
2 = 152 
-*’ = 1 57^ 


<088, 216) 
“<088, 224)“ 
<085, 232) 
<084, 238) 
<080, 246) 
<083, 258) 
-X08.5, 26.8 ).. 




loi P G' 


- Jn-a>3e»3i 



157 RT <D95.. 274) 

158 RT <098. 276) TtA/lAjfO Aa)T SHAO»)w% 


$RLIM RRER 

COLOR: 

GRflV: 

PERCENT: ( 


P AA^A 


0 _ 

" 032 
01 . 02 


_ 1 _ 

064“ 
00. 68 


2 

-096 - 

01. 34 


3 

128' 
03. 74 


4 

■"160^ 
38. 46 


5 

""19 2 ' 
51. 48 


6 

~224~ 
02. 62 


7 

~2“5^" 
00. 54 


$RUN I 

»RUN I 

SCAN 

Z«233 

Z=236 

Z=233 

Z=235 

Z-231 

2 = 232 

Z = 228' 

Z=232 

Z=236 

Z=233 


MRRKER 

HISTOR 


:i50, 

150, 

157, 
159, 
159, 

158, 

151, 
:i49, 
146, 
140, 


276)" 

270) 

266) 

258) 

252) 

242) 

242)' 

254) 

262) 

272) 





b^rS> ■■ I 1 

H iQi2 

lA/i P 2-0 r/A 


Z = 040 
Z = 020 
2 = 022 
2 = 031 
2^013 
Z = 020 
2=028 
Z = 004 
2=^010' 
2 = 013 


RT < 
RT < 
RT < 
RT < 
RT < 
RT < 
RT < 
RT ( 
RT "< 
RT < 


118, 

117, 

118, ’ 
118, 
116, 
114, 
112, 
110 , 
107, 
102 , 


258) 
250) 
240)' 
236) 
228) 
224) 
216) 
232) 
242f 
250) 









iM M r-J tM rsj rsj ru r>i rj im 









APPENDIX G 


MAPPING OF THE 1973 MISSISSIPPI RIVER FLOODS FROM ERTS 
Remote Sensing and Water Resources Management 
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